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We study topological phases of interacting systems in two spatial dimensions in the absence of 
topological order (i.e. with a unique ground state on closed manifolds and no fractional excitations). 
These are the closest interacting analogs of integer Quantum Hall states, topological insulators and 
superconductors. We adapt the well-known Chern-Simons K-matrix description of Quantum Hall 
states to classify such 'integer' topological phases. Our main result is a general formalism that 
incorporates symmetries into the K-matrix description. Remarkably, this simple analysis yields the 
same list of topological phases as a recent group cohomology classification, and in addition provides 
field theories and explicit edge theories for all these phases. The bosonic topological phases, which 
only appear in the presence of interactions and which remain well defined in the presence of disorder 
include (i) bosonic insulators with a Hall conductance quantized to even integers (ii) a bosonic analog 
of quantum spin Hall insulators and (iii) a bosonic analog of a chiral topological superconductor, 
whose K matrix is the Cartan matrix of Lie group Eg. We also discuss interacting fermion systems 
where symmetries are realized in a projective fashion, where we find the present formalism can 
handle a wider range of symmetries than a recent group super-cohomology classification. Lastly we 
construct microscopic models of these phases from coupled one-dimensional systems. 
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I. INTRODUCTION 

The recent discovery of topological insulator j3EH has 
lead to a renewed interest in phases of matter that are not 
described within the usual Landau paradigm of symme- 
try breaking and order parameters— Topological insula- 
tors and topological superconductors, like integer quan- 
tum Hall states, are gapped in the bulk but differ from 
trivial phases in the topology of their electronic states. 
They are characterized by gapless edge excitations that 
reflect the bulk topology. A new aspect of Z 2 spin-orbit 
topological insulators is the role of symmetry (time re- 
versal in this case) in protecting the topological distinc- 
tion. In the absence of symmetry, the topological proper- 
ties, such as gapless edge states, are generally destroyed. 
As with the integer Hall effect, topological insulators 
and superconductors can be described in terms of non- 
interacting particles. A complete classification of such 
free fermion topological phases that are stabl e to d isor- 
der, in all spatial dimensions, has been achieved 18 * 19 !. The 
remaining outstanding questions for fundamental theory 
have to do with interacting systems. 

Interacting topological phases have been studied at 
two levels. First, the stability of the non-interacting 
classification to interactions has been examinecP^HUl l n 
some cases interactions reduce the number of topologi- 
cal phaseiP^H, i.e. two topologically distinct phases of 
free fermions be continuously connected, via intermediate 
states that involve interactions. The second possibility, 
of interactions leading to new phases, not possible within 
non-interacting particles has also been studied. Largely, 
these have attempted to extend the analogy between in- 
teger and fractional quantum Hall states, to topologi- 
cal insulators. Thus, fractional topological insulators 
have been theoretically discussecPSHSS along with lattice 
realizations of fractional quantum Hall insulators^HUl 
These phases are topologically ordered - in that they 
involve fractional excitations in the bulk and feature 
ground state degeneracies that depend on the topology of 
the space on which they are defined. They are also char- 
acterized by a finite topological entanglement entropy 
(TEE) in the ground state. In contrast, integer quantum 
Hall and topological insulators (and superconductors), 
despite being topologically distinct, are not topologically 
ordered. Bulk excitations are essentially like electrons 
or groups of electrons, and the ground state is unique 
when defined on a manifold without boundaries. The 



TEE vanishes for these phases. Henceforth we shall re- 
fer to gapped phases without topological order as being 
short range entangled (SRE) statesPH It seems appro- 
priate to define interacting "integer" topological phases, 
as new topological phases without topological order, but 
which only appear in the presence of interactions. Do 
such phases exist? And if so, how can they be studied? 

In one dimension, topological order is absent, and all 
topological phases found are "integer" (or SRE) phases. 
They include examples like the Haldane (or AKLT) state 
of gapped spin-1 chains^. Using the matrix product rep- 
resentation of gapped statesP™lthey are argued to be 
classified by projective representations of the symmetry 
group (G) or equivalently by the second group cohomol- 
ogy TL 2 (G,C) of symmetry group G. In higher dimen- 
sions, such rigorous results are not available. Neverthe- 
less, new work indicates that these are also amenable to 
theoretical study. Recently, Chen, Gu and WerP^ have 
proposed that higher dimensional group cohomology de- 
scribes d — 2,3 dimensional interacting topological states 
without topological order. For example, bosonic systems 
were studied, where there are no topological phases in 
the absence of interactions. With interactions, topologi- 
cal phases were predicted in two (and three) dimensions, 
without topological order. While Chen et alP^ restrict 
attention to the non-chiral subset of these states (i.e. ones 
that do not have a net imbalance of left and right movers 
at the edge of a two dimensional system) protected by 
symmetry, Kitaev^has also considered chiral states. Ex- 
plicit examples of such phases in special cases has been 
gj ver (48i49j However predictions in the general case rely 
on writing Wcss-Zumino-Wittcn terms for generalized 
sigma models. While this is a powerful approach, the 
physical meaning of the phases that are predicted are 
obscure. For example, the nature of edge excitations in 
these phases is not apparent. Moreover, a knowledge of 
the Borel-group-cohomolgy machinery is required, which 
is mathematically sophisticated even by the standards of 
the field. 

i^-matrix formulation: Here we take a completely 
different and simpler approach to the problem, focusing 
on the case of two spatial dimensions. We rely on the K 
matrix formulation of quantum Hall states, a symmetric 
integer matrix that appears in the Chern-Simons action: 
(h = 1. and summation is implied over repeated indices 
[x,v,\ = 0,1,2): 



4ttS C s = J d 2 xdt Y,i,.j e ^\a I ^]i,jd v a 



(1) 



While this has been utilized to discuss quantum Hall 
states with Abelian topological order, here we show that 
it is also a powerful tool to discuss topological phases 
in the absence of topological order. The latter requires 
|detK| = 1 (i.e. K is a unimodular matrix) . The bulk 
action also determines topological properties of the edge 
states. For example, the signature of the K matrix (num- 
ber of positive minus negative eigenvalues) is the chirality 
of edge states - the imbalance between number of right 
and left moving edge modes. Maximally chiral states 
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1. No Symmetry. 

Z classes. Chiral Edge States. 
Quantized Thermal Hall Conductivity 



/ // / 

2. U(l} Symmetry. 3. Time reversal & U(l) [ha[Ee Symmetry: 

Zclasses. Non-chiral Edge. Z 2 classes. Non-chiral Edge. 

Quantized Hall Conductivity: 

Figure 1: Summary of some simple 'integer' bosonic topo- 
logical phases. 1) A chiral phase of bosons (no symmetry 
required). An integer multiple of eight chiral bosons at the 
edge is needed to evade topological order, leading to a quan- 
tized thermal Hall conductance n xy /T = 8nLo m units of the 

TT 2 k 2 

the universal thermal conductance Lq = 3h B . These are 
bosonic analogs of chiral superconductors. (2) A non-chiral 
phase of bosons protected by U(l) symmetry (eg. charge con- 
servation). Distinct phases can be labeled by the quantized 
Hall conductance a xy — 2n<ro, which are even integer multi- 
ples of the universal conductance ao = q 2 /h for particles with 
charge q. These are bosonic analogs of the integer quantum 
Hall phases. (3) a non-chiral phase stabilized in the presence 
of time reversal and U(l) charge conservation symmetries, 
the same symmetries used to define quantum spin Hall (topo- 
logical) insulators. A Zi topological classification is obtained, 
although bosonic time reversal that squares to +1 is involved. 



have all edge modes moving in the same direction. Phys- 
ically the fluxes e^ LuX d^a v are related to densities and 
currents of bosons of different flavors. 

Strategy and results: Let us briefly review our strat- 
egy and results. Although we mainly focus on non-chiral 
states, we begin by looking for maximally chiral states 
of bosons without topological order. These are bosonic 
analogs of the integer quantum Hall effect or chiral su- 
perconductors of fermions. It is readily shown that the 
smallest dimension of bosonic K matrix that yields a 
maximally chiral state is 8. This is consistent with the 
prediction of Kitaev, derived from topological field the- 
ory. Here we explicitly construct a candidate K matrix 
for this state, corresponding to the Cartan matrix of the 
group Eg. 

We then consider non-chiral states of bosons, with 
equal number of left and right moving edge modes. In 
the absence of symmetry we argue that there are no non- 
trivial topological phases with |detK| = 1. However, 
the presence of a symmetry can lead to new topologi- 
cal phases. The main result of this work is a scheme 
to classify topological phases with | detK| = 1, that are 
protected by symmetry. 

Given a particular symmetry (e.g. time reversal, charge 



conservation etc.) we study distinct ways in which the 
symmetry can act on the elementary quasiparticles. The 
symmetries are realized by a set of symmetry transforma- 
tions on elementary quasiparticles, which form a (faith- 
ful) representation of the symmetry group G. Distinct 
realizations are potentially different phases - analogous 
to the space group classification of crystals. However, an 
additional requirement to realize a non-trivial topological 
phase is the existence of symmetry protected edge states, 
i.e. either the edge is gapless, or if it is gapped, it must 
spontaneously break the symmetry. Note, an internal 
symmetry can never provide such protection to a purely 
one-dimensional system - hence the edge states enjoy 
these special properties by virtue of their connection to 
the bulk topological phase. We will call such phases Sym- 
metry Protected Topological (SPT) phases following the 
terminology of Ref. 47. To access these states, we sup- 
plement the Chern-Simons action with insertion/removal 
of "local" quasiparticles that are bosonic and have trivial 
mutual statistics with any other excitation. Symmetry 
imposes additional, and crucial, constraints on the possi- 
ble terms. The set of these symmetry-allowed perturba- 
tions can be used to analyze if stable edge states exist or 
not. For the most part we restrict our attention to 2 x 2 
K matrices. 

Remarkably, this simple analysis yields the same set of 
interacting topological phases as the group cohomology 
classification of Chen et al for a large set of symmetry 
groups G (see table |H|. For example, bosonic phases with 
G = U(l) are classified by an integer which is just the 
quantized Hall conductance in units of 2Q 2 /h where Q 
is the unit of boson charge. When U(l) is broken to a 
discrete subgroup Z„ , the set of topological phases is also 
reduced to Z„. Similarly both schemes find a Z2 classi- 
fication of bosonic insulators with conserved charge and 
time reversal symmetry (G = U(l) x Zj, the semidirect 
product ensures this is the usual relation between charge 
and time reversal), the analog of fcrmionic quantized spin 
Hall insulators, despite the fact that the time reversal op- 
eration is "bosonic" and squares to +1. An advantage of 
the present formulation is that the edge states of these 
phases are explicit - typically being non-chiral c = 1 con- 
formal field theory (CFT) when gapless. Moreover, be- 
ing cast in the familiar Abelian Chern-Simons form, it 
is amenable to further investigation using standard field 
theory methods. We focus on symmetries (such as time 
reversal) that are realized locally. Spatial symmetries 
such as translation invariance, inversion etc. will be left 
for future work. Since we do not make any assumption 
about spatial uniformity - the topological phases we find 
are well defined in the presence of disorder. 

A disadvantage of our method is that it is less suited 
to discuss non- Abelian Lie group symmetries, and we are 
currently restricted to two spatial dimensions, neither 
of which is a restriction for group cohomology theory^. 
Also, our method does not automatically produce a group 
structure for the set of topological states. On general 
grounds, one expects the set of topological phases pro- 
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Symmetry 


Topological Classification 


Comments 


No symmetry (chiral) 


Z 


_Bg state and derivatives with chiral central charge^ c_ = 8n. 


zl 


Zi 


Time reversal symmetry 


U(l) 


z 


Charge conserved. Quantized Hall conductance a xy = 2nq 2 jh with n £ Z 


17(1) x Zl 


z 2 


Bosonic Quantum Spin Hall with charge U(l) and time reversal T 2 = +1. 


U(l) x 


Zi 


[/(I) spin conservation and time reversal. 




z„ 


1/(1) broken down to a discrete subgroup 


Z n XI Z2 


Z( n ,2) 


(a, 6) = greatest common divisor of a and 6. 


Z n x Z 2 T 


Z (",2) 




U(l) x Z 2 


Z X Z| 





Table I: Topological classification of gapped D=2+l dimensional phases of bosons with short range entanglement (no topological 
order). 



tected by a particular symmetry to form an Abelian 
group, which is automatically satisfied in the group co- 
homology classification and in the classification of free 
fcrmion topological phases. We handle this by defining a 
group structure on pairs of phases described within the 
K matrix formulation. With this refinement the group 
structure of the resulting sets of phases is readily deter- 
mined. 

Topological Phases of Interacting Fermions: We 

extend our discussion to classifying topological phases of 
interacting fermion, in the absence of topological order. 
A key difference from the bosonic case is that symme- 
tries are realized projective ly on the fermion fields. We 
compare our results to a recent super-cohomology clas- 
sification of interacting fermion phases^. In addition 
to the relative simplicity of our method, an advantage 
over super-cohomology classification is that we are able 
to handle Kramers time reversal symmetry T 2 = (— 1) W/ 
(Nf is the total fermion number operator). A disad- 
vantage, shared by the super-cohomology classification, 
is that we are not able to capture chiral or nonchiral 
states with odd numbers of Majorana edge modes. As 
expected, we recover the Z2 classification of time rever- 
sal symmetric quantum spin Hall insulators, from this 
interacting formalism as well. We also compare our re- 
sults with the recent work^HUl on topological phases of 
interacting fermions with Z2 X symmetry. 

Microscopic Quasi-ID Realization: Finally, to 
give a deeper insight into the obtained topological phases 
we utilize a quasi-one-dimensional (coupled wir d 51 ! 52 -*) ap- 
proach to construct a candidate state consistent with the 
edge content that emerges from the classification. The 
K matrix approach naturally suggests such a construc- 
tion. This sheds light on various paradoxical results such 
as the fact that there is a bosonic analog of the quan- 
tum spin Hall state although time reversal acts only on 
bosons with T 2 = 1. 

Some aspects of this work are similar in spirit to a 
number of previous works that have discussed the role of 
symmet ry and stability of edge states in various specific 
contexts ^ 8 | 49 | 53| ^l For example, Ref. 12815^1 discussed the 



stability of edge states in fractional topological insula- 
tors. However, the general machinery presented here to 
generate symmetry protected topological states has not 
previously been discussed. 

II. K MATRIX FORMULATION OF 2 + 1-D 
TOPOLOGICAL PHASES 

It is believecPSHUl x matrix provides a complete 
classification of all Abelian fraction quantum Hall (FQH) 
states in 2 + 1 dimensions. In this section we briefly re- 
view the K matrix formulation of Abelian FQH states. 
We then discuss how it can be applied to study states 
without topological order. In particular we point out 
that in the absence of symmetry, fields that have trivial 
(or bosonic) self and mutual statistics will be 'Higgsed', 
and the stability of the edge is examined in the pres- 
ence of these terms. This is then applied to study chiral 
topological phases of bosons that lack topological order. 

A. A brief review of the K matrix formulation 

The low-energy effective theory of an Abelian quantum 
Hall state is captured by:(h = 1) 

Ccs = ^afajdvai - altf + ■■■ (2) 

(summation over repeated indices is assumed). The 
aj t , I = 1, 2, ■ • • , N are internal gauge fields coupled to 
quasiparticles currents jj, and K is a symmetric matrix 
with integer entries. For states built entirely out of un- 
derlying bosons, the diagonal elements of K are all even 
integers, while for those built from underlying fermions 
(electrons), at least one diagonal entry is an odd integer. 

The topological order is also characterized by the K 
matrix. The ground state degeneracy (GSD) on a torus 
can be calculated by quantizing the Chern-Simons theory 
([2| and is given b}^21 

GSD on a torus = | det K| (3) 
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We will mainly be interested in states without topological 
order i.e. with |detK| = 1. 

Quasiparticles are characterized by integer vector 1, 
and couple minimally to the combination X)j^ a u- The 
self (exchange) statistics 8 of a quasiparticle is obtained 
by integrating out the gauge fields: 

9 = ttFK" 1 ! (4) 

while the mutual (braiding) statistics on taking quasipar- 
ticle li around quasiparticle I2 is: 

012 = 2vrl 1 T K- 1 l 2 (5) 

The effective action describing the gapless edge exci- 
tations of a FQH state characterized by K can also be 
derived^ from gauge invariance of Lagrangian Q on a 
manifold with boundary: 

n I,J 

Here Vj,j is a positive definite constant matrix, that is 
non-universal. However, the commutation relations be- 
tween fields is fixed by the first term that is simply the K 
matrix in the bulk. The number of right movers n + and 
left movers n_ are given by the signature of symmetric 
matrix K, i.e. the matrix K has n + positive eigenvalues 
and n_ negative eigenvalues. 

One important question is: are different FQH states 
characterized by different K matrices fundamentally dif- 
ferent? In other words, can two different K matrices 
represent the same phase? This means that two FQH 
states have exactly the same set of quasiparticles but 
these qu asiparti cles are labeled in two different ways. It 
turns ouPE2Ill tnat 

a generic change of label (or change 
of basis) for the same set of quasiparticles is realized by 
the following GL(N,Z) transformation: 

<^E^-><< w eGL(iV,Z). (7) 
j 

Here GL(N, Z) denotes all N x N integer matrix with 
determinant ±1. After this relabeling of quasiparticles 
the K matrix and currents i], transform as 

K -> W T KW, (8) 

Any two K matrices related by such a GL(N,Z) trans- 
formation represent the same state (in the absence of any 
global symmetry). It's straightforward to see that physi- 
cal properties such as the determinant and the signature 
of a if matrix is invariant under such a GL(N, Z) trans- 
formation. 

When there is a U(l) symmetry associated with charge 
conservation, one couples an external U(l) gauge field 
to the conserved U(l) current with charge q via an inte- 
ger vector t = (*!,••■ , £at) t called the charge vectoi^. 



This is incorporated by adding the following term to the 
Lagrangian ^ above: 27r£ cha rgc = -q^vxtiA^duO^. By 
integrating out internal gauge fields {a^} one can obtains 
the quantized Hall conductance: 

a xy = |Vk-H (9) 

Z7T 

and the U(l) charge of a quasiparticle with integer vector 
1 is given by Q = {jt T K _1 l. 

The many-body wavefunction of a multi-layer FQH 
state described by effective theory Q is given b} * 63 * 66 ! 

* K = n (^ ) -4 7) ) Kj ' j ^ e - k(i,|2/4 - (io) 

i<j,I,J 

in a disc geometry. Here zf^ = + \y\ de- 

notes the two-dimensional coordinates of the i-th par- 
ticle in the 7-th layer. Multi-particle pseudopotentials 
can be constructecP^HSSl as ideal Hamiltonians, whose 
zero-energy ground states are the above multi-layer FQH 
states. 



B. K matrix + Higgs formulation 

The Lagrangian ^ seems to have U(1) N symmetry (or 
N conserved currents) due to the existence of N inter- 
nal gauge fields {a^}. When these correspond to bosonic 
excitations (featured by trivial self and mutual statis- 
tics with other quasiparticles) , and in the absence of any 
symmetry, one generically does not expect them to be 
conserved. This can be implemented by adding terms 
to the action Q that create and destroy these bosonic 
particles, which we (in the absence of a better phrase) 
call Higgs terms. To be precise, denote the annihilation 
operator for a quasiparticle of 7-th type as bj and the 
associated creation operator as bj =6}. If an integer 
vector 1 = (?i, • • • , Im) T characterizes a boson, then we 
demand: 

tt1 t K" 1 1 = mod27r (11) 

and 

27rl T KT 1 l' = mod 2tt (12) 

for trivial mutual (braiding) statistics with all other 
quasiparticles 1'. Then, in the absence of symmetry we 
can add a Higgs term 

SC CS = Ci I]/ b 1 / + h.c, Ci = const. (13) 

to the Lagrangian Ccs that condenses the boson with 
vector 1. Note, since this particle has trivial statistics 
we can dispense with the gauge field in this expression, 
whose only role here is to keep track of statistics. One can 
show that the topological properties of the corresponding 
state, such as the quasiparticle statistics Q and ground 
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state degeneracy ^ are not influenced by these Higgs 
terms. Therefore the generic Lagrangian describing a 
2 + 1-D gapped Abelian phase is the following 



C 



2 + 1 



c 



cs 



E (cill# 



h.c. 



(14) 



{1— bosonic} 



Taken at face value, the Chern-Simons theory which at- 
taches flux to particles would require monopole terms to 
account for a change in particle number. We have argued 
above why this may be un-necessarjS2J In any event, as 
we will see below, the only action of the Higgs terms we 
will need is their effect on the boundary - which does not 
suffer from these problems. 

Stability of Edge States: Now, the action of edge 



excitations corresponding to bulk Lagrangian (14 1 is 



9 1 = V 

edge ^^{1— bosonic} 



Sedge — S° dge + ^edgei ( 15 ) 

C\ J dtdxcos(J2ih4>i + a i)- 



where S® d is given in (|6|. The bare action S® dge indi- 
cates the following Kac-Moody algebra: 

[dMx), d y cbj(y)] = 27riKJ 1 J d x S(x - y) (16) 



Notice that each allowed Higgs term ( 13 1 in the bulk has 



a one-to-one correspondence with those on the edge in 



S edge> Le - 



CiY[b 1 / + h.c 



-+ 



Cicos(^ Ijcbj + on) 



(17) 



While all these perturbations are present at the edge a 
more stringent requirement needs to be placed if they 
are to gap out the edge modes. For example, we expect 
a maximally chiral edge, where all modes move in the 
same direction, to be stable even in the absence of any 
symmetry. The required condition can be deduced by 
studying the commutation relation implied by the Kac- 
Moody algebra above for the field 1 T • <fi = J2i^i ( t ) i : 

[1 T • d x cb(x), l T ■ d y ct>{y)] = 2 7 ri(l T K- 1 l)a E <5( a; - y) (18) 

thus, in order to be able to localize this field at a classical 
value, and gap out an edge mode, we require that the 
commutator vanishes i.e. 



i J rn = o 



(19) 



For a maximally chiral state where K 1 is a positive def- 
inite matrix, no non-vanishing vector satisfies this condi- 
tion. Hence the edge states cannot be gapped. Similarly, 
when there are an imbalanced number of right and left 
moving modes, n + ^ n_ the system has a net number of 
chiral modes and we call it a 2 + 1-D chiral phase. In a 
2 + 1-D chiral phase even in the absence of any symmetry 
there will be gapless edge excitations 53 . 

To completely gap out an edge one requires equal 
number of counter-propagating modes, i.e. a non-chiral 



edge. Then, dimension of K matrix N is even and 
N/2 = n+ = n_. Let us call i.e. cos(li T + a\ ± ) and 
cos(l2 4> + a h) independent Higgs terms if and only if: 



h T K-\ 



1 2 T K 1 1 2 



h T K-% 







(20) 



In this case they form a pair of commuting variables ac- 
cording to the Kac-Moody algebra. According to Heisen- 
berg's uncertainty principle, these mutually commutat- 
ing fields {l n T </>|n = 1, 2, • • • } can be pinned at certain 
classical values simultaneously, and consequently their 
associated edge excitations will be gapped. Then, to 
completely gap out the edge one needs a set of N/2 in- 
dependent Higgs terms that are pair wise commuting. In 
the absence of any symmetry, this is typically possible. 
However, in the next section we will see that symmetry 
can forbid some Higgs terms leading to SPT phases with 
non-trivial edge structure. Now let us first consider a 
chiral state of bosons without topological order. 



A chiral bosonic phase without topological 
order: The E% state 



A phase without topological order is characterized by 
a symmetric K matrix with | detK| = 1. A chiral state 
in 2 + 1-D requires the signature (n + , n_ ) of its K ma- 
trix to satisfy that n + ^ n_. Such a state has gap- 
less edge excitations and a non-zero quantized thermal 
Hall conductance^. There are many such examples for 
a fermionic system: e.g. an integer quantum Hall state 
whose K matrix is the unit matrix of size N . On the 
other hand, in a bosonic system without topological or- 
der, the existence of such states is less obvious. 

We therefore seek a K matrix with the following prop- 
erties (i) |detK| = 1 (ii) the diagonal elements K] j 
are all even integers and (iii) a maximally chiral phases, 
where all the edge states propagate in a single direction. 
Then, all eigenvalues of K must have the same sign (say 
positive), so K is a positive definite symmetric unimod- 
ular matrix. 

It is helpful to map the problem of finding such a K 
to the following crystallographic problem. Diagonalizing 
K and multiplying each normalized eigenvector by the 
square root of its eigenvalue one obtains a set of prim- 
itive lattice vectors ej such that Kjj = ej ■ e,/. The 
inner product of a pair of vectors liei and Z^e/ are given 
by I'jKijIj, while the volume of the unit cell is given by 

1/2 

[DetK] ' . The latter can be seen by writing the compo- 
nents of the vectors cts a, square matrix: - [c/]a- 
Then Detfc is the volume of the unit cell. However, 
Kjj = E a k aik a j = (k T k)„. Thus DetK = [Detk] 2 . 

Thus, for a phase without topological order, we re- 
quire the volume of the lattice unit cell to be unity 
[Detfc] = 1 (unimodular lattice). Furthermore, for a 
bosonic state, we need that all lattice vectors have even 
length IjKjjIj = even integer, since the K matrix has 
even diagonal entries (even lattice) . It is known that the 
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minimum dimension this can occur in is eig htP. In fact, 
the root lattice of the exceptional Lie group Ee is the 
smallest dimensional unimodular, even lattice^. Such 
lattices only occur in dimensions that multiple of 



A specific form of the K matrix is: 



/ 2 -1 



K = 









V o 



-1 2 
-1 








-1 



-1 2 






-1 



-1 2 
-1 






-1 



2 

-1 







1 

-1 
2 

-1 




o \ 


-1 




-1 




2 / 



(21) 



This matrix has unit determinant and all eigenvalues are 
positive. It defines a topological phase of bosons without 
topological order, with eight chiral bosons at the edge. 
Note K" 1 can be related to K by a GL(8,Z) transfor- 
mation S T K~ 1 S if we take S — K, and so they are 
physically identical. Thus, even without computing the 
inverse, all particles are seen to have trivial statistics 
{irl' T K~ 1 l = 2-Krn). This K is the Cartan matrix for 
E$, hence the name of the state. 

This result was previously pointed out by 
utilizing the fact that the central charge of the edge states 
(c_ = c — c) of a chiral topological phase are determined 
by the statistics of emergent excitations only modulo 8. 
Thus, phases with trivial statistics are allowed whenever 
c_ = mod 8. 

Combining these phases leads to an integer classifi- 
cation of chiral topological states, which are character- 
ized by a quantized thermal hall conductivity 70 which 
are integer multiples of the universal thermal conductiv- 
ity: lim T ^ ^ = 8^". 



III. INCORPORATING SYMMETRIES IN K 
MATRIX FORMULATION 



In this section we will be interested in incorporating 
global symmetries into the K matrix + Higgs formula- 
tion. This will lead to new symmetry protected topologi- 
cal phases. We only consider internal symmetries, spatial 
symmetries like inversion, translation etc. will not be dis- 
cussed. 

Now let us restrict ourselves to 2 + 1-D non-chiral 
phases with equal numbers of counter propagating modes 
n + = n_ for signature (n+, n_) of matrix K. In the ab- 
sence of symmetry any edge Higgs term that satisfies ( |13[ ) 
can be added. In such a phase there will be no gapless 
edge excitations in the absence of any symmetry, i.e. all 
edge modes will be gapped out by the Higgs terms S], dge . 

However this is not true any more when there are sym- 
metries in the system. In the presence of symmetry, only 



those bosonic quasiparticles which transform trivially un- 
der the symmetry operation can condense. This means 
certain Higgs terms which transform nontrivially under 
the symmetry operation are not allowed and cannot be 
added to effective theory (14) or ( |i5| . 

How do quasiparticles transform under a symmetry op- 
eration? In general the Lagrangian Q and ([Mf should be 
invariant under the symmetry transformation on quasi- 
particle currents {jj }. Notice that when a K matrix is 
acted on by a GL(N,Z) transformation (|8j, it describes 
the same physical state. Only the labels of different 
quasiparticles are changed. Given a state described by a 
certain K matrix, the allowed GL(N,Z) transformations 
W which transform the quasiparticles under symmetry 
must leave the K matrix invariant i.e. K = W T KW. 
Besides any global U(l) phase transformation on the 
quasiparticle annihilation operator bi — > e^^'bj also 
keeps the Lagrangian ^ invariant. Notice that such a 
phase shift S(f>i is defined modulo 2ir due to the quanti- 
zation of quasiparticle number in a compact theory 

Therefore a generic realization of unitary symmetry g 
on the {<fii} fields at the edge (we will not require the 
transformation law in the bulk in the following) is given 

by 



(22) 



where 6(f>j € [0, 2ir) are constants and matrix W 9 £ 
GL(N,Z) satisfies 



K = (W 9 ) T KW 9 . 



(23) 



For an anti-unitary symmetry h (such as time reversal 
symmetry Zj), in general it is realized in the following 
way: 



(24) 



where drf) 1 } E [0, 2ir) are constants and matrix W h € 
GL(N,Z) satisfies 



K = -(W h ) T KW h . 



(25) 



The anti-unitary symmetry operation h is realized by 
the above transformations followed by complex conjuga- 
tion C. Notice that K matrix changes sign under the 
above GL(N,Z) transformation since a Chern-Simons 
term e tll/X aj l d l/ a , l always changes sign under time rever- 
sal. 

Group compatibility conditions: It may appear 
we have wide latitude in determining how the generators 
of a symmetry group act on the quasiparticles in our 
theory. However, there is an important constraint. For 
an arbitrary symmetry group G, the multiplication rule 
of its group elements is completely determined by certain 
algebraic relations of the group generators {g\, Q2, ■ • ■ }■ 



■4{«a> = n 



9 a 



(26) 
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where e is the identity element of group G and {n a } are 
all integers. A bosonic quasiparticle (which satisfied Eqn. 
Ill is a physical excitation and must transform trivially 
under the identity element. Thus all boson insertion op- 
erators satisfying ( 11 ) and (l2"| should be invariant under 
the symmetry operation A{ na y. 

-4{n.} : J^ihfa ^J2i l i<Pi mod 2?r, (27) 
V 1 satisfying fK" 1 ] = mod 2. 

These algebraic requirements serve as constraints to the 
possible GL(N,'Z) transformations W 9a and U(l) phase 
rotations {S(/)j a } and we shall call them group compati- 
bility conditions for the edge states described by effective 
theory (15 1 and their associated bulk topological phases. 



In the case of bosonic phases without topological order, 
K is a symmetric unimodular matrix whose diagonal el- 
ements are all even integers. Then, any integer vector 1 
satisfies the conditions (111 and (12 1, i.e. all quasiparti- 
cles of a bosonic SRE phase are bosons. Therefore the 



group compatibility conditions ( 27 1 for symmetry trans- 
formations are simplified as 

Under A{ Ua } = e : (f>j — »• 0/ mod 2n, (28) 
1 = 1,2,..- ,N. 

By solving these algebraic equations we can find 
out all sets of inequivalent symmetry transformations 
{W 9a , 8(f> 9a } for generators {g a } of group G. 

As an aside we note that for phases with topological 
order, or in the presence of fermionic quasiparticles, sym- 
metries are realized projectively^H Then, even the iden- 



tity elements ( 26 1 can induce a nontrivial transformation 



on quasiparticles. 

Gauge Equivalence: A question naturally arises: do 
different symmetry transformations represent different 
SPT phases? We answer this question in two parts. First 
we comment on the equivalency or inequivalency between 
two sets of symmetry transformations {W 9a , S(j) 9a }. No- 
tice that one can always change the label of quasi- 
particles by a GL(N,Z) transformation X (as long as 
X T K.X = K), or perform a global U(l) gauge transfor- 
mation 4>i — > tfii+A(f)i which keeps Lagrangian |2| invari- 
ant. Under such a "gauge" transformation the symmetry 
operations {W 9a , <50f a } will transform as 



W 9 -» X~ l W 9 X, 



(29) 



x- 



if X £ GL{N, Z), X T KX = K. 



where r) = ±1 if g is a unitary (anti-unitary) symmetry. 
If two sets of symmetry operations {W 9a , 5(j)j a } associ- 
ated with group G generated by {g a } are related by the 
above gauge transformation ( |29| , then these two sets of 
symmetry operations are essentially identical. 

Edge Stability and Criteria for SPT Phases: 
With a set of symmetry transformations {IU 9a , Scf) 9a } one 



allowed in the presence of the symmetry group G = {g}. 
In general only those Higgs terms transform trivially un- 
der the symmetry transformations {W 9a ,6(f> 9a '} are al- 
lowed, which induce certain allowed set of edge pertur- 
bations. To determine the fate of possible gapless modes 
at the edge, one considers terms that commute with each 
other and with themselves, i.e. terms involving variables 
^2jh<j>i that satisfy conditions (19) and (20 1. Then one 



can simultaneously minimize these terms like classical 
variables. Note, scaling dimensions of these edge terms 
are immaterial to this discussion. 

Criterion for Trivial Phase: If there is a set of inde- 
pendent Higgs terms {C a cos(\J (p + a a )} allowed by sym- 
metry, so that any other variables l T cj> on the edge is ei- 
ther a linear combination of these bosonic variables {1^ </>} 
or doesn't commute with every condensed bosonic quas- 
particles in {1^0}, then the edge of the system will be 
completely gapped in the presence of independent Higgs 
terms {C a cos(rj0 + a a )}- When the independent Higgs 
terms {C a cos(1^0 + a a )} are simultaneously minimized 
on the edge, the elementary bosonic variables {v^0} 



v a = (v a i, ■ ■ ■ ,v a N ) 



In 



Z,cd(l a ,l,l a ,2, ■ ■ ■ >L,n) 



, V a(30) 



can determine what Higgs terms in (|14| and ( 15 ) that are 



will all condense and be localized at various classical val- 
ues (Va4>) = B a (gcd is short for greatest common di- 
visor). Notice that if the set of independent elementary 
bosonic variables {vj^} are invariant under any symme- 
try transformation i.e. 

Under VgeG: £> a ,/<M "> ( 31 ) 



then the edge states can be all gapped out without break- 
ing the symmetry G at all. We call such a non-chiral 
SRE phase a trivial phase since in general it doesn't sup- 
port gapless edge states in the presence of symmetry G. 
These principles will be illustrated by examples in detail 
in the following section. In comparison to the aforemen- 
tioned trivial phase, a nontrivial SPT phase has a gapless 
edge structure which cannot be gapped without breaking 
the symmetry G. The only two possible situations for 
the edge structure of such a phase are: 

(i) Gapless: The maximal set of independent Higgs 
terms {C a cos(l^</> + a a )} allowed by symmetry cannot 
gap out all the edge states. In other words, on the 
edge there exist at least one variable lip that commutes 
with all the condensed bosonic variables {1^0}. Hence 
this degree of freedom \<j) remains gapless even in the 
presence of all the symmetry-allowed independent Higgs 
terms {C a cos(l^(/) + a a )}. An example is the bosonic 
SPT phases protected by Z7(l) symmetry in 2 + 1-D, as 
will be discussed in detail later. 

OR 

(ii) Spontaneous Symmetry Breaking on the Edge: Al- 
though all the edge states will be gapped in the presence 
of independent Higgs terms {C a cos(1^0+o; a )} allowed by 
symmetry, not all of the associated elementary bosonic 
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variables {vj0} in (30) are invariant under symmetry 



transformations {r) g W g , 8cf> 9 }, Vg £ G where G is the 
symmetry group. This me ans at least one elementary 
bosonic variable v„ (j> in ( 30 1 would transform nontrivially 



under symmetry group G. Therefore in order to gap out 
the edge by condensing all the independent elementary 
bosons, one has to spontaneously break the symmetry 
on the edge. An example is the bosonic SPT phase pro- 
tected by U(l) x Z2 symmetry (or by Z2 symmetry) in 
2 + 1-D, as will be shown later. 

Group structure of phases protected by sym- 
metry group G: In general, the set of different phases 
that appear in a topological classification are expected 
to form an A belian group, as proved for non-interacting 
fermion d 19 ! 74 ! and conjectured for interacting bosonic 
systems^ (since group cohomology classification leads 
to an Abelian group). How does this group struc- 
ture appear within our K matrix formulation? Let 
{V G [K, {W 9a ,Scf> 9 <-}]} to denote the set of trivial and 
SPT phases in the presence of symmetry G. Here 
^g[K, {W 9a , S(f> 9a }} represents a phase with matrix K 
in action (15 1 and transformation rules {W 9a ,84> 9a } for 
symmetry group G = {g a }- We would like to attach to 
this set a group product. 

A natural Abelian product rule © of two phases 
$ G [K,{W 3a ,5<l) a «}] and * G [K, {W 9 « ,S~(/) 3a }} is to take 
their matrix direct sum: 

* G [K, {W 9 ° , 8^ }] © * G [K, {W 9 « , <f/ a }] 

= * G [K © K, {W 9 * © W B ',S^' © 6~<fi 9a }}. (32) 

This seems to suggest that one cannot obtain a full clas- 
sification of all different SPT phases when restricted to 
a K matrix with a fixed dimension. However notice that 
two phases \I/ G and ^> G can be identified as the same one 
if © eg = i.e. adding a trivial phase (denoted 
by e G ) to ^f G yields the phase We can use this fact 
to reduce the dimensions of K matrix by throwing away 
the "trivial" parts of the edge structure, which can be 
gapped without breaking any symmetry. 

The identity element in the group {^g} corresponds 
to the trivial phase e G = ^ G [K Q ,W 9a = lN X N,5(/) 9a = 
0], where Ko can be any N x N unimodular symmetric 
matrix corresponding to a non-chiral SRE phase in 2 + 1- 
D. It edge states can be gapped out without breaking the 
symmetry of group G. 

We can also define the "inverse" of a phase 
* G [K, {W 9a , 8(f) 9 -}] in the group to be 

* G [K,{W B %8<i> 9 *}]- 1 = y G [-K,{W B %8<i> a «}}. (33) 

i.e. by changing the sign of its K matrix we obtain the 
inverse of a phase. This is simply because we can always 
gap out the edge of 

^gIK^W^,^}]- 1 ®y G [-K 7 {W 9 «,8<j) 9 «}} = e G 

without breaking the symmetry. Consider the variables 
{</>/} of phase [K, {W 9a , 8<j) 9a }] and {0/} of phase 



* G [-K, {W 9a , 8<j) 9a }] on the edge. Then edge pertur- 
bations such as: {C a cos [l„ (0 — <f>) + a a ]} will not be 
affected by the phase factors 8(f> 9a . We can then con- 
dense a set of independent elementary bosonic variables 
{l^((j>— <j>)} satisfying (19 1 and (20). This can be readily 
shown when either W = I (arbitrary K) or when K is a 
2x2 matrix (arbitrary W). This is simply because {cj)i} 
and {<fri} satisfies the Kac-Moody algebra with opposite 
K matrices. No symmetry will be broken by condensing 
these bosons, when a proper set of vectors {l a } are cho- 
sen. To check if two putative SPT states are the same 
phase or are different phases, we use the above group mul- 
tiplication rules to combine one state with the inverse of 
the other and check if a trivial phase results. If so, these 
two are the same SPT phase. 

Now that the identity element, the inverse of an el- 
ement and the multiplication rules are defined, we can 
identify the group structure. We will perform this anal- 
ysis below to clarify the connection between phases gen- 
erated by our formalism. 

Miscellaneous Considerations: We focus on the 
cases where {W 9a ,84> 9a } form a faithful representation 
of symmetry group G, where G is the symmetry group 
of the system's Hamiltonian. The case when it forms 
an unfaithful representation, i.e. when more than one 
group element acts like the identity on all quasiparticles 
(this set of group elements form an invariant subgroup 
H), actually corresponds to the faithful representation 
of the quotient group (G/H). Thus, studying faithful 
representations suffices. This is discussed in more detail 
in Appendix [D] In some cases we will find that a solution 
for symmetry transformation may not be realizable in a 
theory with local fields (for example one that exchanges 
a pair of fields that are canonical conjugate of each other 
in a c = 1 edge), which will then not be included in the 
minimal set of SPT phases. For phases that are reported 
in Table |TT] we have checked that they have symmetries 
that can be realized starting from microscopic degrees of 
freedom (such as from the coupled wire construction). 

In the following we illustrate the above principles, first 
by classifying bosonic SPT phases in 2+ 1-D, whose topol- 
ogy is protected by a certain symmetry group G. 



IV. AT-MATRIX CLASSIFICATION OF 
BOSONIC SPT PHASES 

In this section we will focus on bosonic non-chiral SRE 
states in the presence of certain symmetry, i.e. bosonic 
non-chiral SPT phases in 2 + 1-D. They are described 
by a symmetric, unimodular K matrix, whose diagonal 
elements are even integers and with the same number of 
positive eigenvalues and negative eigenvalues (n+ = n_). 
Therefore the dimension of matrix K must be even and: 

dctK = (_i) dim ( K )/ 2 for a non-chiral SRE phase. 

In this section we restrict ourselves to a 2 x 2 K matrix. 
SPT phases that are necessarily described by K matrices 
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of a larger size, may be missed by this restriction. How- 
ever our results are internally consistent and also capture 
at least all the topological states of the group cohomol- 
ogy classification of bosonic SPT phases in Ref. H7J It 
appears that a K matrix of dimension 2 is sufficient to 
represent and classify SPT phases in 2 + 1-D in many 
cases. The reason behind this unexpected success is the 
following: although we focus on SPT phases described 
by a K matrix of size 2x2, when analyzing the group 
structure formed by SPT phases with symmetry G we 
need to multiply two phases together by a direct sum of 
their K matrices. Consequently we are in fact consider- 
ing K matrices of size 2n x 2n obtained from direct sums 
of original 2 x 2 K matrices. Therefore it's actually not 
surprising that many of the bosonic SPT phases in 2 + 1- 
D can be described and classified by a 2 x 2 K matrix 
and associated symmetry transformations. 

As proved in Appendix [S] a 2 x 2 K matrix with de- 
terminant — 1 for a bosonic system (see n = 1 in theorem 



where e is the identity operation. Time reversal symmc- 



( B3 1) is always equivalent to the standard form a x — 
10 I 



1 0; 



by certain GL(N, Z) transformations (a a , a 



x,y,z are Pauli matrices). In the following we always 
choose the 2x2 matrix K = a x to represent a generic 
bosonic non-chiral SRE state. In the following we use 
general principles discussed earlier in section (IIP to 



study non-chiral bosonic SPT phases with different sym- 
metries. Note that the GL(2,Z) transformations X that 
keeps K = a x invariant under ^ are X = ±^2x2, ±o" x . 
For such a non-chiral bosonic SRE phase, the unper- 
turbed edge theory is: 

^edge = ^J dtd * (Mldxfa + dtfadsfa (34) 

- YjI,j Vi,jd x (f>id x <pj) 

where Vi j is a positive definite constant matrix, as dis- 
cussed in section [II B| This implies the following commu- 
tation relations (Kac-Moody algebra) for the edge fields 

i>i,<M : 



[d x <f>i(x),d y <f>i(y)] = [d x <p 2 (x),d y <p 2 (y)] = 0, 
[d x <pi(x), d y fa{y)] = 2rtid x 5(x - y). 



(35) 



In the absence of any symmetry, the edge states in ( 34 ) 



can always be gapped, since a set of independent Higgs 
terms satisfying ( |T9| and ( 20 1 on the edge can be choose 
as either {Ci cos(l<fii+ai), I £ Z} or {Ci cos(l<f>2+on), I £ 
Z}. These added terms destroy the edge states. All de- 
grees of freedom on the edge are gapped when variable 
cj>\ (or (f>2) is localized at a classical value. 
Now, let us consider various symmetries. 



A. Z% symmetry: Zi class 



symmetry (time reversal) is generated by T. The 
algebra ( 26 ) which defines symmetry group is 



try T is implemented (following rules ( 24 1 ) by a matrix 
W T £ GL(2, Z) and a vector of U(l) phase changes S(j) T 



(defined modulo 2tt) satisfying the constraint ( 28 1 and 
( 25 1 for a bosonic SRE system 



x2- 



T\T 



KW T = K; 



(37) 



Ej Wf,jWj = mod 2tt, 7=1,2. (38) 



The only GL(2,Z) matrix solutions to (37 1 are W T = 
±<7 2 . Notice that W T = a z and W T = — a z are related 
a gauge transformation (29 1 X = a x . Therefore one can 
always choose 



W 



by a proper gauge fixing. Then the constraint (38 1 be- 
6 T = mod 2tt and it leads to 



comes (hx2 - <?z 



b![ = H.27T mod 2ir, n,2 = 0, 1. 



Under a gauge transformation A(j)j in ( 29 1 the compact 
U(l) phase shift 5(f> T transforms to 8<t> T — \I2-x2 +er z )A</>. 
As a result we can always choose a gauge so that 

ScpJ = mod 2tt. 

So a generic bosonic non-chiral SPT phase in the presence 
of time reversal symmetry T has symmetry transforma- 
tion 

{W T ,5^ T } = {aJ ° )} n 2 = 0,l. (39) 
\n2-Kj 

Since each bulk Higgs term has a one-to-one correspon- 
dence with that on the edge, hereafter we'll only write 
down those Higgs terms C\ cos(^ J li4>i + ct\) on the edge. 
In the case of symmetry the allowed Higgs terms are 

S ldge = 'El 1 >Q,l 2 Clfdxdt COS^i + hfc + <*{) 

+ cos(-^i0i + ; 2 02 + nihn + a{) . 

For ri2 = 0, 1 the allowed Higgs terms are different, 
e.g. cos(02) terms are allowed for ri2 — but not allowed 
for n,2 = 1. Thus there is a distinction between these 
states. However, this is not a topological distinction as 
argued below. For both 712 = and n 2 = 1 cases we 
can write the same set of symmetry-allowed independent 
Higgs terms formed by mutually commuting operators 
(|35): i.e. 



S ldge = T,h C h,0 J dxdtCOS (h<j)l(x,t) 



(40) 



(36) 



Thus, in both cases the edge can be gapped, therefore 
they belong to the same trivial phase. If the variable <\>\ is 
localized at expectation value eg. (4>\) — by the Higgs 
terms, all excitations on the edge would be gapped but 
the time reversal symmetry (4>i — ¥ — </>i) is not broken by 
this expectation value. 
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B. U(l) symmetry: Z classes 

The elements of U(l) group can be labeled as U$ where 
9 £ [0, 2tt) and the identity element is Uo. The multipli- 
cation rule is given by 



Ug 1 Ug 2 — U(e 1+ g 2 



od 27t) 



(41) 



and therefore A generic form of symmetr y t ransforma- 
tions {W 9a ,54> 9 j a } satisfying constraint (28 1 for U(\) 
group is 



{W u > = Iax2, 



et 




1 }, {tiM) = l- 



where (£i , £2) denotes the greatest common divisor of in- 
tegers t\ and £2- Notice that only when (£1,^2) = 1 
the above symmetry transformations {W Ul> , <j) Ue } form 
a faithful representation of symmetry group U(l), which 
is what we assume hereP^ 



Here t' = Kt = 



is nothing but the charge vector 



defined in the context of K matrix formulation of a FQH 
state (see section |II A I . As proved in Appendix [C] for a 
bosonic SRE phase with K = o~ x , an arbitrary charge 
vector t with (£i,£2) = 1 is equivalent to the standard 
form 



, qez. 



by certain GL(N, Z) gauge transformations. Therefore 
the inequivalent symmetry transformations under con- 
straint ( 28 1 for a bosonic non-chiral SRE phase with U(l) 
symmetry are 




{W u » = I 



2x2; 




}, qez. 



The associated symmetry-allowed Higgs terms are 



el 

edge 



J2iez^i J dxdt cos 



(42) 



Due to Kac-Moody algebra ( |35| , they do not form a set of 
independent Higgs terms satisfying ( 19 1-( 20 ) unless q = 



0! When q — the Higgs term can localize the variables 
02 to expectation value (02) = const., which gaps the 
edge excitations without breaking f/(l) symmetry (0i — ¥ 

01+0). 

Now we'll determine the group structure formed by 
these phases. Let's label a phase by 



[q] = * u , 1) \K = a tB ,{W u ' =/ 2x2 , 



kUe - 




variables {0'i,02}- When they are put together the fol- 
lowing independent edge terms are symmetry allowed: 
J2iezCi cos [l(4>i — (/)[)+ a 1] and their associated elemen- 
tary bosonic variable is 0i — (f>[ 
dent Higgs terms Yliez ^1 cos \} 
allowed by U(l) symmetry. Therefore the edge states 
will be fully gapped without breaking the G = U(l) 
symmetry, by condensing elementary independent bosons 
{0i — (/)[, 4>2 + 4>2}- Therefore we have 



If (72 = —qii indepen- 
2 + 02 ) + a 'i] are a ^ so 



V q e Z. 



(43) 



On the other hand if q± q-2. the bosonic variable = 
02+02} cannot be gapped since cos |j(02+02)+ a !] terms 
are not allowed by symmetry. Now the new variables 
describing the gapless edge structure can be chosen as 
{01 = 0i, 02 = 02 + 02 } satisfying Kac-Moody algebra 
(35 1. Notice that under U(l) they transform as W Ue — 
I 2X 2 and 



Al + 52, 



Therefore we have the multiplication rule of the group 
formed by phases [q] with U(l) symmetry: 




[<7i] © [92] = [qi + 92] 



(44) 



This means phases [q] labeled by different integer q's are 
different phases in the presence of U(l) symmetry, and 
they form nothing but the integer group Z! Any phase [q] 
with q =/= corresponds to a nontrivial SPT phase, whose 
gapless edge states cannot be gapped without breaking 
U(l) symmetry. 

There is a simple physical reason underlying these ob- 
servations. The Hall conductance is a physical invariant 
that distinguishes these different phases: 



tVK 



-1+' 



K 



2q. (45) 



C. U(l) X Z2 symmetry: Z2 classes 

In the presence of both "charge" U(l) (group elements 
Ug) and time reversal Z 2 symmetry (generator T), the 
extra algebraic relation in addition to (36 1 and (41 1 is 
given by 

U-gT = TUg 

since the charge U(l) symmetry doesn't commute with 
time reversal symmetry. The algebraic relations for 
U(l) xi Zl are 



TUgTUg = e. 



(46) 



where e-un) — [0] is the trivial phase. Consider two in addition to (41). The corresponding constraints (28 1 
states [q±] with edge variables {0i, 02} and [52] with edge for symmetry transformations {H /T , 5<fi T } and {W Va = 



12 



9t} are 



{h*2~W T )(5<f + 0t) = 



mod 2tt, V 9(47) 



and ( 37|38 I . Just like in the case of Z^ symmetry, again 



using gauge transfor mat ion X = o T i n ( 29 1 on e can fix 
W T = o z to satisfy (37). Solving (38) and (47 1 we have 
and 64 



t-2 = and 64>i, = nn, n = 0, 1. Using gauge transforma- 
tion A4> in (29 1 one can fix 84>\ = 0. Hence the inequiva- 
lent symmetry transformations for 17(1) X Zj symmetry 
is 



W Ue = I- 



2x2, 




(48) 



W =a z , 64> T = n = 0, 1. (49) 

\ nir 



The symmetry-allowed independent Higgs terms are 



ci 

^ edge 



lei 



Ci / dxdt 



cos(Z0 2 ) + cos(Z02 + nln) 



(50) 



Apparently they consist of only the 4>2 variable and 
hence are independent of each other. When n = all 
cos(/02)j I € Z terms are allowed in (50) and it corre- 
sponds to the trivial phase. Variable 4i can be local- 
ized at value (</> 2 ) — and the edge excitations will be 
gapped without breaking any symmetry. On the other 
hand when n = 1 only cos(Z0 2 ), I =even are allowed in 
(50 1 and it corresponds to the nontrivial SPT phase. If 



4>2 is localized at any value the edge becomes gapped. 
However time reversal symmetry T will be broken since 
under T we have 2 — ^ 4>2 + n - 

Now let's analyze the group structure formed by these 
phases. We denote the two phases with K = a x and 
symmetry transformations (49 1 where n — 0, 1 as [0] and 
[1]. [0] = e v ^^ Z T is the trivial phase. Consider two 
copies of nontrivial SPT phases [1] put together: they 
have edge variables {^1,^2} and {^,^> 2 }- Apparently 
one can always condensed the pair of elementary inde- 
pendent bosons {4>i — 4>\i §2 + ^2} on the edge, and the 
edge states will be fully gapped without breaking any 
symmetry. Therefore we have 



[i]e[i] = [o]. 



(51) 



Clearly [0] and [1] form a Z 2 group. As a result n = 
and n = 1 label the Z 2 classes of bosonic non-chiral SPT 
phases for U(l) x symmetry. 

Note, the generator of charge U(l) Ug = e mB does not 
commute with time reversal which involves complex con- 
jugation which sends i — > —i in the exponential. Thus 
U-qT = TUg, which implies time reversal and charge 
conjugation are combined via the semi-direct product 
U(l) xi Zj. However, if the U(l) was associated with spin 



rotation about S z for example, of an integer spin system, 
the relation with time reversal would be that of a direct 
product [7(1) x Z2, since now UgT = TUg. This com- 
pletely changes the topological classification and leads to 
no nontrivial phases as shown in Appendis jE 1| 



D. Zn symmetry: Zjv classes 

Denoting the generator of Zn group as g, the algebraic 



structure ( 26 1 of Z^ group is given by 



9 N = e. 



(52) 



The corresponding constraints ( 28 ) for symmetry trans- 
formations {W 9 , 64 9 } are 



(W 9 r=I 2x2 , (W^ a x WB = a. 



mod 2-7T. 



(53) 
(54) 



In the following we discuss the cases of N being an odd 
and even integer respectively. 



1. N — odd integer: Zjv classes 

It's straightforward to check that the only solution 
7 2 x2- So the solutions to (54) have 



to (531 is W 9 



the general form of 64 9 — where t = ( 1 ) and 

w 

(h,t 2 ) = (k,N) = 1 for k,h,t 2 = 0, l,-- - ,N- 1. Here 
we require = (k,N) — 1 so that the transfor- 

mations ( |55[ ) form a faithful representation of symme- 
try group G — Z^. Making use of theorem (CI I, we 



can always reduce an arbitrary "charge vector" t with 



(*1,*2 



1 to its standard form 



and hence the 



J 



inequivalent symmetry transformations {W 9 , 64 9 } for 
Zn, N = odd symmetry are 



W 9 = hx2, 
{k,N) = l, 




0,1,---,N-1. 



(55) 



It's easy to show that the Higgs terms allowed by sym- 
metry don't depend on k and they are 



ci 

edge 



+ql 2 =0 mod N 



Ci J dxdtcos(l T + aiX56) 



Notice that a Higgs term labeled by vector 1 is allowed 
only if l\ + ql 2 = mod N. Apparently when q = this 
is a trivial phase with a set of independent Higgs terms 
being X^ez Q J dxdt cos(Z0 2 + a;), and the variable </) 2 
can be localized at any value without breaking the Zn 
symmetry. Since for different k values in transformations 
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( 55 I the symmetry- allowed Higgs terms are exactly the 



same, we believe different k values correspond to the same 



phase and we will assume a representative k = 1 in ( 55 1 
hereafter. 

To analyze the group structure of these states, let's de- 
note the various phases with symmetry transformations 
( 55 I under Zjv symmetry by 



[q} = ^ ZN [a x ,{W B = I 2x2 , 



27T 1 

N [n 



}}■ (57) 



where [0] = [N] = &z N is the trivial phase. Again con- 
sider two states [qi] with edge variables {(j>i, <fi 2 } and [q 2 ] 
with edge variables {(/)[, 2 }. Completely in parallel with 
the discussions for U(l) symmetry, it's straightforward 
to show that 



[N-q], [qi] [q 2 ] = [<7i + q 2 mod N]. (58) 



Therefore different phases [q] with q = 0,1, ■ ■ ■ , N — 1 
form a Zjv group. There are Z^v classes of different phases 
labeled by q — 0, 1, • • • , N — 1 in the presence of Zjq 
symmetry, when N — odd. 



2. N = even integer: Zn classes 



W B = 



Now the inequivalent solutions to ( 53 1 

±-'2x2 j ±0"a!- 

(i) For W 9 = I 2x2 we have exactly the same solutions 
as ( p)5| ) in N = odd case, and hence Z^r different classes 
of bosonic non-chiral SPT phases. All these Z^ phases 
can be realized by coupled wire construction, as will be 
discussed in section IVII 

(ii) For W 9 = —I 2x2 one can always choose a gauge 

so that 8<fi 9 — 0. The symmetry- allowed Higgs terms are 
J2y !COs(l T 0) and this describes nothing but the trivial 



phase as [q = 0] in (55). Its edge states can be gapped 
out without breaking any symmetry. 



3. N = even integer: other solutions 



(iii) For W 9 — a x , the gauge inequivalent solutions 



to (54 1 are 
We 



i> 9 = ^jr 



1 



with k = 0, 1, • • • , N - 1. 



require (k, N/2) = 1 so that these transformations 
{W 9 , 5<p 9 } form a faithful representation of symmetry 
group Zn, N =even. The symmetry-allowed Higgs terms 
are 



S ldge = E 2 fe(il+i 2 )=0 mod JV C l I dxdt ' ( 59 ) 

cos(l T (/> + ai) + cosihfo + l 2 ^ x + ?^±M + ai) 



One can verify that they all corresponds to nontrivial 
SPT phases, whose edge states cannot be gapped without 
breaking the symmetry. More precisely, variables <pi and 
</>2 cannot be localized simultaneously since they do not 
commute, and if only one variable (say (f>i) is localized 
the symmetry g will be broken since 4>i 4>2 + ^w- 
Their symmetry transformations are summarized as 



W 9 = a x , Scj) s 




(fc,f) = l. (60) 



But it's not clear how to realize these phases in a mi- 
croscopic model or what group structure they form. We 
label these phases as [a x , k] = ^ z N [o'x, {W 9 = a x , 5<f> 9 = 

^ (])}}■ 



1 



(iv) For W 9 



ful" solutions to (54) are 5<f) 9 = 



with k = 



a x , the gauge inequivalent "faith- 

0, 1, • • • ,N -I ancTl$, N/2) = 1. The symmetry-allowed 
Higgs terms are 



2nk 
N 



Cl 

edge 



E 2k (h-l 2 )=0 modN C lJ dxdt - ( 61 ) 

cos(l T + ai) + cos(-ii^2 - hfii + 27rfc ^~' 2 ^ + oci) 



If we label these phases as [— a x ,k] = ^>z N [o'x, {W 9 



| }], it's easy to see that 



We discuss below additional representations of the 
symmetry group that appear for this particular case, 
which we believe are unphysical for a SRE phase with 
no ground state degeneracy on a toruJ^. These re- 
quire interchanging the two edge fields 4>i,4> 2 , which have 
very different character when they describe fundamental 
bosons (one is like the phase field, and is compact, while 
the other is related to the integrated density) . Therefore 
we believe it is unphysical to exchange them. Also, unlike 
for the other symmetry transformations, a microscopic 
model with this realization of symmetries was not found. 
Finally these additional phases are not naturally accom- 
modated into a group structure. These points taken to- 
gether lead us to drop them from the final list of topo- 
logical phases with this symmetry. 



[<y x ,k] 



-fx, k] 



(62) 



This is because once a [a x , k] state with edge variable 
{<j>i, ^2} and a \—<j x , k] state with edge variable {4>i, 4> 2 } 
are put together, one can always condense bosons 
{</>i + <pi, 4> 2 — </> 2 }, and the edge will be gapped without 
breaking Zn symmetry. 

To summarize, no matter N = odd or N — even, there 
are Zjv different bosonic non-chiral phases in the presence 
of Zn symmetry. They are characterized by different 
symmetry operations (55 1 associated with Zn generator 
g and symmetry-allowed Higgs terms (56 1. All these Z^v 



phases can be realized in coupled wire construction as 
will be shown in section IVT1 
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Besides, when TV = even there are extra solutions (60) 



to constraint ( 54 I for symmetry transformations associ- 
ated with Zn symmetry. However the physical realiza- 
tion of these states and their group structure are not 
clear. 



E. Zn x Z 2 symmetry 



The generators of Zn x Z 2 symmetry group are g for 
Zn and T for Z 2 satisfying the following algebra 



,7V 



TgTg = e. 



(63) 



The associated constraints on symmetry operations are 
W 9 W T W 9 W T = J 2x2 , (64) 
(I 2 x2-W B W T )(5(f> B + W B 5cj> T ) = (°) mod27r. 



in addition to ( |37p8| and (|53|54l 



1. N = odd integer: Zi class 



The gauge inequivalent solutions to these constraint 



equations are ( 39 I and 




(k,N) = l. 



(65) 



Let's label these phases by [712, k] where n 2 is defined in 
( 39 1 . Notice that when n 2 = in ( 39 I one can always 



destroy the gapless edge excitations by localizing variable 
02 without breaking any symmetry (under T we have 
4>2 4>2 + ti 2 7t). Similarly when k = in ( |65| the edge 
can be gapped out by localizing bosonic variable 0i. So 
n 2 = or k = both correspond to the trivial phase. 
On the other hand when n% — 1, the symmetry-allowed 
Higgs terms are 

= Ei Ci / dsdif cob(JVIi^ + / 2 2 + ai) (66) 



- C0S(— Nli<f)i + / 2 02 + Ctl + i 2 ^2 7r ) 



At first sight it seems the edge states cannot be gapped 
without breaking the symmetry, i.e. neither 0i nor 2 
can be localized due to symmetry. However when a state 
[n 2 = l,k =/= 0] with edge variable {0i,0 2 } is put to- 
gether with a trivial state [1,0] = e ZN>>Z T with edge 

variable {0'i,0 2 }, the edge can be fully gapped by con- 
densing bosons {N<f>i + (/)[, cf> 2 — N(f>' 2 } without breaking 
the Z N symmetry (N =odd). Therefore [l,fc] © [1,0] = 
e z N »z T an d are a U trivial phases, which in gen- 

eral doesn't have gapless edge states. As a result there 
is no nontrivial SPT phases in the presence of symmetry 
Z N x Z%, N =odd. 



2. N — even integer: 

When N — even we always have W 9 — ±/ 2 x2- 
(i) For W° — I 2X 2, the gauge inequivalent faithful so- 
lutions to the constraint equations are ( 39 I and 



(k,N/2)=l, n = 0,l. (67) 



Let's label various phases with symmetry transforma- 
tio ns (|39| and (67) as [k, ri2, n] where 712 = 0, 1 is defined 
in (|39[). When k = variable <j>\ can be localized without 
breaking any symmetry and it is the trivial SPT phase. 
When rt 2 = n = the variable 4>2 can be localized and it 
is the trivial phase again. Therefore 



[0,n 2 ,n] = [A, 0,0]. 



for N = even. In the following we analyze the group 
structure formed by these states. 

Following discussions in section |III| we can obtain the 
inverse of a phase by merely changing the sign of its K 
matrix. Now let's put together a state [k,ri2, n] with edge 
variable {</>i, ^2} is put together with a state [k', n' 2 , n'}^ 1 
with edge variable {$1, $2} , we can condense the follow- 
ing independent bosonic variables {k'4>\— k<j>'i, kfc— k' fy} 
and destroy the gapless edge states if (k,k') = 1. The 
associated Higgs terms will not break the Zn x Z 2 sym- 
metry if k'n—kn' = mod 2 and fcn 2 — k'n'o = mod 2. 
As a result [k, n 2 , n] © [A;', n 2 , n']" 1 = e,, 



'2 ^ 



For (k,k') — 1 : [fe, n 2 ,n] = [fc', ; 



''2) ' 



if fe'n — kn' = mod 2, kri2 — k'n' 2 = mod 2. 

Therefore we have [2fc + l, n 2 ,n] = [1, n 2 , n]. For fc =even 
on the other hand, we know that N/2 must be odd 
since (k,N/2) — 1 for a faithful representation. Then 
we can choose k' = and condense independent bosons 



{^4>i — 4>i, 4>2 — t° destroy all edge state. No sym- 

metry will be broken by doing so. Hence we showed 
[2fc,7i2,n] = [0,n2, n] = &z N -»zJ- Consequently the 
only three nontrivial SPT phases are [1, 1, 0], [1, 0, 1] and 
[1,1,1]- 

Similarly by putting together a state [1, 712,71] with 
edge variable {(^1,(^2} is put together with a state 
[1,77,2,77,'] with edge variable {<^,<^ 2 }, we can a lways lo- 
calize bosonic variable <pi — and gap out part of the 
edge. What is left on the edge are described by variables 
= 0i, 02 = 02 + 02 }• They obey Kac-Moody algebra 
(35 1 and transform as a [1, n 2 + n' 2 ,n + n'] state. Hence 



we've shown that 



[1, n 2 , n] © [1, n' 2 , n'] = [1, n 2 + n' 2 , n + 



(69) 



Since 77, n, 2 = 0, 1 are both Z2 integers, so clearly all 
different 4 states [1, n 2 , n] form a group. Consequently 
there are 3 nontrivial SPT phases labeled by k = 1 and 
[n 2 ,n] = [0,1], [1,0] or [1,1] in @ and ((67 K 
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(ii) For W° 



(29 1 so that 8(f) 9 = 



12x2 we can always choose a gauge in 
V 



Higgs terms added is 



,0, 



From constraint equations one 



can derive W T = a z and (72x2 ± W T )8cj) T 
we have 



Til, 712 = 0, 1. 




and 



(70) 



Apparently if m — variable <f>i can be localized without 
breaking any symmetry, and similarly if Ti2 = variable 
<f>2 can be localized without breaking the symmetry. For 
the nontrivial SPT phase with m =712 = 1, the edge 
states cannot be destroyed without breaking any symme- 
try. If we label this SPT phase by [ni = 1, 7i2 = 1], one 
can show that the group structure formed by states with 
symmetry transformations W 9 = —I is the integer group 
Z i.e. {[1,1]", n € Z}. However the above symmetry 

transformations {W 9 — — ^2x2, <5</> 9 = ^j^J } do not 

correspond to a faithful representation of Zn X group 
for N =even unless N = 2. And it is not clear whether 
the states with symmetry transformations W 9 = —I can 
be realized in a physical bosonic system. Therefore we 
won't include the states with symmetry transformations 
W 9 = — I in the minimal set of topological phases with 
Zjv x Z2 symmetry. 

To summarize, there are 13^ classes of different non- 
chiral bosonic SRE phases (including one trivial phase 
and three SPT phases) with Zjv x Zj symmetry when 
TV = even. To compare, the classification and analysis 
of bosonic SPT phases with Zjy x Zj symmetry (the 
direct product of Zjv and Zj in contrast to the semi- 
direct product discussed here) is shown in Appendix |E 2| 



A-MATRIX CLASSIFICATION OF 
FERMIONIC SPT PHASES 



According to theorem ( B3 1 in Appendix |B|a 2 x 2 K 
matrix with determinant —1 for a fermionic system is 

always equivalent to the standard form ^ M ~ cr z 

by certain GL(2,Z) transformations. In the following 
we always choose the 2x2 matrix K = u z to repre- 
sent a generic fermionic non-chiral SRE state. In the 
following we use general principles discussed in section 
|III| to study non-chiral fermionic SPT phases with differ- 
ent symmetries. Note that the only GL(2,Z) transfor- 
mations X that keeps K = a z invariant under (|8| are 
X = ±/ 2 x2j ifj. For such a non-chiral fermionic SRE 
phase, its "bare" Chcrn-Simons effective theory with no 



£k = ^\ald v a\ - ald v a\) - ]T atf (71) 



4tt 

in the bulk and 



1=1 



Sedge = h J dMx {MMl ~ d^Afc (72) 

-J2i,j v i,J d x(f>idx(f>j) 

on the edge where V/,j is a positive definite constant ma- 
trix, as discussed in section [Tl B| The Kac-Moody algebra 
satisfied by fields {4>i, 4>2} on the edge is 

[d x <f>i(x),d y <t>i{y)] = -{d x 4> 2 (x),d y (f> 2 (y)] = 
2irid x 5(x-y), [d x <j>i(x), 9 y 2 (y)] = 0. (73) 

In the absence of any symmetry, a set of independent 
Higgs term satisfying (19) and (20 1 on the edge can 
be chosen as either {C; cos(Z</>i + l(f>2 + I € Z} or 
{CiCOs(l<fii — l<t>2 + Q-i), I € Z}. All degrees of free- 
dom on the edge will be gapped once "bosonic" variable 
01 + 4*2 ( or 4>2 — 4*1 ) is localized at a classical value by the 
Higgs terms. Again the two bosonic variables (pi+(j>2 and 
— 4>2 cannot be localized simultaneously, according to 
Heisenberg uncertainty relation implied by Kac-Moody 
algebra (73 1. 



There is an intrinsic difference between fermions and 
bosons: i.e. only bosonic quasiparticles can "condense" in 
a bosonic / fermionic system described by a local Hamilto- 
nian. This means in a bosonic system, all quasiparticles 
are bosons and should transform trivially under identity 
element e = FT g" a of group G as shown in ( 28 ) . In a 



fermionic system on the other hand, any bosonic quasi- 
particle consists of an even number of fermions and is 
always invariant if every fermion creation (annihilation) 
operator obtains a minus sign. Therefore under symme- 
try transformation {W e , 5(fi e } corresponding to identity 
element e of the same symmetry group G in any local 
fermionic system, only those Higgs terms satisfying con- 
dition (11|-(12| should transform trivially. This can be 



generalized to a more universal situation where anyonic 
quasiparticles are present (|detK| > 1): under identity 
element of group G only loc a l op erators (i.e. Higgs terms 
cos(l T + ai) satisfying ( 11 )-( 12 1 which condense bosonic 
quasiparticles) should transform trivially. This means 
with the same symmetry group G — {g a }, the symme- 
try transformations {W^ a ,<5<^f a } of a bosonic SRE state 
form a faithful representation of group G, while symme- 
try transformations {W^ a , <5</>^ a } of a fermionic SRE state 
(or more generally a gapped Abelian phase containing 
fermionic and anyonic quasiparticles) form a projective 
representatiorP^ of group G = {g a }- And for these sys- 
tems the identity element e in group compatibility condi- 
tions (26 1 and (27 1 doesn't always correspond to a trivial 



transformation on the fermionic (anyonic) quasiparticles. 

In a fermionic non-chiral SRE phase with K = a z here, 
it is easy to verify that such a bosonic quasiparticle is 
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Symmetry 


Minimal Topological Classification 


Comments 


Zl (no symmetry) 


Zi 


No symmetry (fermion parity always conserved) non-chiral phase 


Zl x Z{ 


Zi 


Time reversal symmetric superconductor 


17(1) x Zl, T 2 = 1 


Z 2 


Bosonic quantum spin Hall insulator of Cooper pairs 


17(1) x Zl, T 2 = P f 


zi 


Including fermionic quantum spin Hall insulator 


(7(1) x Zl x Z| 




?7(1) spin conservation and time reversal. 


Z 2 x Z| 


z 4 


Superconductor with Ising-type symmetry 


Z 4 


z 2 


bosonic ^-symmetric SPT phase of Cooper pairs 


Z 2 x x Z| 


Z 4 x Z 2 


discussed in Ref. 21 23 with Zs classification 


Z 4 x 


zi 


Superconductor with Z4 spin symmetry 


Z 4 X Z^ 


zi 


Time-reversal-symmetric charge-4e superconductor 



Table II: Topological classification of gapped D=2+l dimensional non-chiral phases of fermions with short range entanglement 
(no topological order). Here Pf denotes fermion parity, which is always conserved. Note, states with an odd number of 
right(left)-moving Majorana edge modes (such as class Dili topological superconductor) are not captured in this formalism. 



labeled by any vector 1 



h 



satisfying 



mod 2 



(74) 



Locality requires that only Higgs terms cos(l T (f> + a\) sat- 



isfying (74 1 can be added to bare action ( |71| and ( |72| , i.e. 
fermions are not allowed to condense. Hence in the ab- 
sence of any symmetry, all these Higgs terms cos(l T cf>+a\) 



satisfying (74) are allowed and should be added to a 
fermionic non-chiral SRE state. As a result, the iden- 
tity element e (no symmetry) in a fermionic system is 
implemented by the following generic form of symmetry 
transformations 



'2x2, 



J7/7T 



(75) 



where r\j — 0,1. Notice that the above symmetry 
transformations are invariant under any gauge transfor- 
mation (29). When rjt = 1 the fermionic operators 
a = 1, 2 obtains a minus sign, corresponding 



exp[it 



to the fermion number parity operation Pf = (—1 
When rjf ~ on the other hand, every fermion remains 
invariant and it corresponds to the actual identity ele- 
ment e/ of the symmetry group Gf (including Pf) for 
the underlying fermions. If we incorporate the fermion 
number parity P f into the symmetry group G / , one eas- 
ily notices that z£ = {ef,Pf} is always a normal sub- 
group of fermion symmetry group Gf, which means Pf 
is involutory (P f 2 = e/) and central in Gj^. 

If we always incorporate fermion number parity Pf 
into fermion symmetry group Gf, then a fermionic sys- 
tem with symmetry Gf is naturally related to a bosonic 
system with symmetry G = Gf/Z^ (since is a nor- 
mal subgroup of G f the quotient group Gf/Z^ can be 
defined) . Physically this means if fermions with symme- 
try Gf pair up to form Cooper pairs (which are bosons), 
these bosonic Cooper pairs have symmetry G — GfjZ 2 . 



Different SPT phases of bosonic Cooper pairs (where 
fermions are confined) with symmetry G — Gf/z£ are 
necessarily different fermionic SRE phases with symme- 
try Gf. Hence the different classes of fermionic SRE 
phases with symmetry Gf must contain all different 
bosonic SRE phases with symmetry G = Gf/Z^ MS H 
subseiPS. To be more precise, the group Hf,(G f /Z£) 
formed by different bosonic (non-chiral) SRE phases with 
symmetry G = Gf/z[ is always a subgroup of Hf(Gf), 
the group formed by different fermionic (non-chiral) SRE 
phases with symmetry Gf. 

Before discussing specific examples of non-chiral topo- 
logical phases, we point out that SRE chiral phases of 
fermions are readily obtained (e.g. integer quantum Hall 
states) so we skip their discussion. It should also be noted 
that our formalism currently is restricted to topological 
phases in which the gapless edge, when present, has in- 
teger central charge (e.g. c = 1 in many cases). So chi- 
ral and non-chiral Majorana modes (e.g. of a p x + ip y 
superconducto r, or the D=2+l class Dili topological 
superconductor^^, with a pair of counter-propagating 
Majorana modes) are not captured in our formulation. 



A. Gf/Z( = {e} => Gf = Z( symmetry: Zi class 



If we choose G 



Gf/Zl 



{e} is the trivial group, 



then the fermion symmetry group is Gf = Z^ . The gen- 
erator of Z<± i.e. fermion number parity operator is 



(-1)"/ 



(76) 



where Nf denotes the total fermion number. The ex- 
istence of Z 2 symmetry is a basic requirement for any 
fermionic system described by a local Hamiltonian. Sim- 
ply speaking P f guarantees that one single fermion can- 
not condense like the bosons: only a bosonic conglomer- 
ate containing an even number of fermions can condense 
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and obtain a non-vanishing expectation value. A general 
form for such a bosonic quasiparticle in a fcrmionic sys- 
tem is labeled by a integer vector 1 satisfying condition 

(fill. 



This means is more like a constraint for fermionic 
system due to locality, rather than a true "symmetry". 
As discussed earlier it is implemented by nothing but the 
non-trivial [nf = 1) realization for e in G = Gf/Z^: 



'2x2, 



Sep 1 



(77) 



It guarantees that in the absence of any symmetry, all 
Higgs terms cos(l T + a\) satisfying (74 1 can be added 
to a fermionic non-chiral SRE state and these are the only 
terms that can be added. Notice that (77 1 is invariant 



under any gauge transformations ( 29 I for fermions with 

2 



K = er 2 . Apparently we have P 



f 



if 



i.e. 



number parity acting twice would yield the identity op- 
eration for fermions. 

The Z 2 symmetry-allowed Higgs terms are all terms 



associated with bosonic quasiparticles ( 74 I 



ci 

edge 



£{l 1= i a mod 2} Ci J dxdt cos(F0 + 0-0.(78) 



Apparently the bosonic variable fa + fa (or fa — fa) 
can be localized at a classical value and the edge will 
be gapped. Physically fa + fa corresponds to the pair- 
ing between right mover and left mover, while fa — fa is 
backscattering between right and left movers. They are 
both allowed in the absence of any symmetry. This de- 
scribes the (trivial) Zi class of non-chiral fermionic phase 
with Z 2 symmetry. 

Since fermion number parity Pf is always realized by 
(751, in the following we'll not specifically mention this 



symmetry but only requires Z 2 to be a normal subgroup 
of the full symmetry group Gf of fermions. And we use e 
to denote the identity element in the "bosonic" symmetry 
group G — Gf /z£ ■ Therefore in the fermion system e 
can be either e f (all fermion operators keep invariant) or 
Pf (all fermion operators change sign). 



on symmetry transformations {W T ,S(f) T } are 

(W T ) 2 = 7 2x2 , (W T ) T KW T = -K; (80) 

(l2x2-W T )5cj) T =nTv(l\ mod2vr, (81) 

mod 2tt. (82) 



1 



where ^TiVTPf — 0, 1. Notice that fermion number par- 
ity symmetry Pf is always implemented by ( 77 1 , inde- 
pendent of the gauge choice. Here with K = a z the gauge 
inequivalent solutions to (80) is W T = a x . Then solv- 
ing (|81 h and (|82|) we get t]t = r\TP, and 5<j) T = [ ^ | . 

U T 7T/ 



Therefore the ineq uiva lent symmetry transformations for 



the fermion Z 2 x z l group is (177 1 



W 



Vt = 0, 1. 



(83) 



And the symmetry-allowed Higgs terms are 

E{/ 1=i2 mod 2} °i I dxdt cos(l T + ai ) 



ci 

edge 



■ cos(-l T fi x + l 2 rjT^ + Qi^ 



(84) 



where t)t = 0, 1. It turns out the above Higgs term 
always describes the same trivial SPT phase no matter 
i]t = or 1: e.g. variable fa + fa can be always local- 
ized at an expectation value (fa + fa) = 772^ /2 by the 
Higgs terms, and the gapless edge states will be destroyed 
without breaking any symmetry. Therefore ( [72] ) together 
with (84 1 describes the (trivial) 1\ class of fermionic non- 



chiral SPT with x Z% symmetry, no matter T 2 

(Vt = VTP f = 0) or T 2 = P f (r? T = Vtp s = 1)- Note, 
the Z 2 classification of free fermions in class Dili with 
these symmetries is missed by this classification. The 
reason of course is that we are not able to describe Ma- 
jorana modes (the class Dili topological superconductor 
has a counter-propagating pair of Majorana modes with 
central charge c = 1/2) within the if-matrix formulation. 



Gf/Zl 



Gf = Z2 x Zl symmetry: Zi class 



C. G f /Zl = (7(1) xi Zl symmetry 



In the presence of time reversal symmetry with gen- 
erator T, the algebraic structure of full symmetry group 

G f/ Z 2 = z 2 is gi ven b Y 



T = e. 



(79) 



Here in our notation e can be either e f , the identity ele- 
ment for fermions or Pf. the fermion number parity op- 



eration. From (27) this leads to the following constraint 



By labeling the U(l) group elements as Ug, the alge- 



braic structure of G = U(l) x Z£ group is given by 



T = TUgTUe = U le=t 



(9=0 mod 2tt) 



= e. 



(85) 



in addition to (41 ). Again here in our notation e can be 



either identity element e f for fermions or fermion number 
parity Pf. Again a general form of symmetry transfor- 
mation for Ug is given by 



b u » = et 



(86) 
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and the we have the following constraints for the symme- 
try transformations 



2nt = ??[/(i)7r I J mod 2n, 



(h x2 - W T )(6t + 5(t) T ) = rjn [ \ mod 2tt. 



(87) 



in addition to (80]81 ). The last line in (85) are automat- 
ically satisfied. The gauge inequivalent solutions to these 
constraint equations are 



t = (t+5^) ' ' 



W T = <r x , 



t G Z; (89) 
mod 2?r. (90) 



where rju(r\, rj = 0,1. Notice that the fermion par- 
ity generator Pf is always a subgroup of t/(l) since 
f + i)[/(i)/2^0. Hence Z| is always a subgroup of £7(1) 
group associated with fermion number conservation. If 
7] = 1 we have T 2 — Pf, while r\ — corresponds to 
T 2 = 1. 



L G/ = (7(1) xi ™£/i T 2 = 1: Z 2 c/asses 



When rj = the algebra of symmetry group 
[7(1) xi Zf is 



T[7 e T[/ e = [7 



(6=0 mod 2tt) 



(91) 



where e / is the identity element of symmetry group G t 
for fcrmions. The symmetry-allowed Higgs terms associ- 
ated with symmetry transformations ( 86 I and ( 89 ) are 



J2iez.Cl J dxdt cos l(fa - fa) + m 



for 77 = 0. Hence in r\ = case there is only one trivial 
phase, since independent bosonic variable fa — fa can be 
localized at a classical value by the Higgs terms and the 
edge will be gapped without breaking any symmetry (un- 
der T we have fa — fa — > <fii — fa — rjir). Meanwhile notice 
that for a bosonic system with G = [7(1) X Zj symmetry 
(T 2 =l) there are Z 2 classes of different phases. Hence the 
nontrivial bosonic SPT phase of Cooper pairs (fcrmions 
are confined) protected by [7(1) x Zj symmetry form a 
nontrivial SPT phase of fermions with Gf = [7(1) xi Zj 
symmetry. As a result there are Z 2 classes of differ- 
ent fermionic (non-chiral) SRE phases in the presence of 
[7(1) X Z-f symmetry with T 2 = 1. The Z 2 classification 
comes purely from the bosonic SPT phases (bosonic QSH 
insulator) of Cooper pairs in the molecule limit where 
fcrmions are confined. 



2. G f = U{1) xi Zl with T 2 = P f : I? 2 classes 



When rj = 1 the algebra of symmetry group Gf — 
U{1) x is 

T 2 = TUgTUg — [7(0 = o m od2,)=Pf- (92) 

And the symmetry-allowed Higgs terms on the edge are 



S ld Q e = Eiez °i I dxdt cos 



2/( 



o 7 



for 77 = 1. This corresponds to the nontrivial SPT phase, 
whose edge cannot be gapped without breaking the sym- 
metry. We use [rj\ with rj — 1 and rj = 1 to label 
these two phases. Now let's examine the group struc- 
ture {^u(l} >(l z T xz , } f° rme d by these states. 

The trivial state is labeled by identity element 
e u(i)>iz T - If we P 11 ^ two [?y = 1] states together, we 
can gap out the edge states without breaking the sym- 
metry so [1] © [1] = Bxj{X)-aZ^' They form a Z 2 group. 
As a result, [77 = 1] and [1] © [1] label the Z 2 classes 
of fermionic non-chiral SRE phases in the presence of 
[7(1) xi Z% symmetry with T 2 = Pf. When [7(1) sym- 
metry corresponds to fermion charge conservation, these 
two different phases are nothing but the trivial band in- 
sulator and Z2 topological band insul ator (qu antum spin 
Hall insulator) of fermions in 2 + l-l j 54 l 75 l 76 l. Again the 
bosonic SPT phases of Cooper pairs (where fermions are 
confined) with G — [7(1) x symmetry gives rise to 
another Z 2 classification. Hence in total there are Z| 
classes of different SPT phases with [7(1) x symme- 
try (T 2 = P f ). 



D. G f /Zl = 7/(1) x Zl symmetry: Zi class 

The algebraic structure of G = G f /z£ = [7(1) x Zj 
group is given by 



T 2 = TU-gTUg = U { g =0 



mod 2-tt) 



(93) 
(94) 



in addition to (41). The associated constraints (|27j) for 
symmetry transformations (86 1 and {W T ,S<j> T } are (80 
81L (|87| and 



(hx2 - W T )S0 T + 0(I 2x2 + W T )t = 7]w ( J j mod 2tt. 



Again the last line of (93) is automatically satisfied. The 



gauge inequivalent solutions to these constraint equations 
are 



t = (i+^P) ' ' 



-1 



W T = a x , 



t G Z; 



mod 2ir. 



7 7 7r 



(95) 



(96) 
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The associated symmetry-allowed Higgs terms are 

S ldge = E;ez °i I dxdt cos [Kfii + <h) + on 
for 77 = and 



S l edge = Eiezl dxdt ( c i sin [( 2/ + + 0a)] 
+£>, cos [2l(0i + fa)] 

for 77 = 1. In both cases variable </>i + 02 can be localized 
at a classical value by the Higgs term, so the gapless edge 
states will be destroyed without breaking any symmetry 
(under T we have fa + fa — > r/Tr — fa — fa). So there is a 
(trivial) Zi class of fermionic non-chiral SRE phase in the 
presence of U(l) X x Z;{ symmetry. Note that T 2 = ej 
if 77 = and T 2 = P t if 77 = 1 . Since there is no bosonic 
SPT phases with U(l) X Z 2 = Gf/Z^ symmetry, there 
are no new fermionic SPT phases coming from bosonic 
SPT phases of Cooper pairs. 



E. Gf/Z{ = Z2 symmetry 

The generator g of Z 2 symmetry satisfies the following 
algebra 



g = e. 



(97) 



Here e stands for either identity element e f for fcrmions 
or fermion number parity P t . This algebraic constraints 
(27 1 for symmetry transformations {VK 9 , 5 fa 1 } are 



(W 9 ) 2 = I 2x2 

(hx2 + W 9 )Sfa = V n 



(W 9 ) T K.W 9 = K, (98) 

V 



mod 2-7T. 



(99) 



where r\ = 0, 1. The gauge inequivalent solutions of (98) 
are W 9 = ±/2 X 2, ±cr z . In the following we analyze those 
cases with W 9 — ±^2x2 and the discussions about cases 
with W 9 — ±(T Z is put in the Appendix. It's not clear 
to us now whether the transformation laws with W 9 — 
±a z can be realized in a microscopic model, therefore we 
didn't include these cases in the minimal set of different 
fermion SRE phases with G/ symmetry. 

(i) For W 9 = — ^2x2 the gauge inequivalent solution 

-ft: 



to (991 is 



and 77 = 0. A set of indepen- 



dent symmetry- allowed Higgs terms can be either (102) 
or (103) with ai = 0, V / G Z. Hence it corresponds 



to the trivial phase, whose edge can be gapped without 
breaking any symmetry. 



(ii) For W° = ^2x2 the inequivalent solutions to (99) 



are 




*i,*a,»? = 0,1. (100) 



The symmetry allowed Higgs terms are those cos(l T <; 
a\) terms satisfying 



l\t\ + /2^2 



T) = mod 2 



(101) 



and the condition (74) for local operators 
forward to verify that when t\ + t% 



It's straight- 
i] = mod 2 a set 



of independent symmetry-allowed Higgs terms satisfying 
fll9l is 



Sldge = Q J dxdt cos [l(fa + fa) + a,] (102) 

and the edge states can be gapped without breaking the 
Z 2 symmetry (under g we have fa + fa — > fa + fa + 
(<i + t 2 + ?7)7r), if variable fa + fa is localized at any 
classical value. Similarly when t\ — t 2 = mod 2 a set 
of independent symmetry-allowed Higgs terms satisfying 
(fl9) is 



el 

edge 



Eiez D i I dxdt cos [l(fa - fa) + Pi] (103) 



and the edge states will be gapped without breaking the 
Z 2 symmetry (under g we have fa — fa — > fa — fa + 
(ti — t 2 )ir). if variable fa — fa is localized at any value. 
They all correspond to the trivial phase. Notice that 
when 77 = we have g 2 — ej while g 2 = Pf if 7; = 1. 
This corresponds to the following two different symmetry 
groups: 



1. Gf = Z2 x Z 2 symmetry: Z4 classes 
Intrinsic fermion phases 

This means 77 = and ti—t 2 = \ mod 2. The algebra 
of symmetry group Gf is 

g 2 = e f . (104) 

When [77, 1 1 , £2] = [0,0,1] or [0,1,0] a set of indepen- 
dent symmetry-allowed Higgs terms satisfying (19) can 
be chosen to be either 



°edge 



I dxdt cos [ 2Z (^i + 02)+ ai] (105) 



= J2 Di [ dxdt cos - + a] ( 106 ) 



They correspond to two nontrivial SPT phases, where 
the edge cannot be gapped without spontaneously break- 
ing the Z 2 symmetry. Let's label these two states as 
[77, ti,t 2 ] = [0, 0, 1] and [0, 1, 0]. Notice that when we put 
a [0,0, 1] edge with variables {fa, fa} together with a a 
[0,1,0] edge with variables {fa, fa}, the edge can always 
be gapped by condensing e.g. independent bosonic vari- 
ables {fa + fa, fa — fa} and no symmetry will be broken. 
Hence [0, 0, 1] © [0, 1, 0] = e. z ^ z s and we have [0, 0, 1] = 
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[0, 1, O]^ 1 . On the other hand, if we put four [0, 1, 0] states 
with edge variables {4> R , <fi%, a = 1, 2,3,4} together, the 
edge can be gapped without breaking the symmetry, by 
localizing the following independent bosonic variables: 



+ 



and hence 



As a results we have [0,1,0] = e z z f 

[0,1,0] 3 = [0,0,1]. Therefore all different fermionic 
phases form a Z4 group. 

To summarize, with Z 2 symmetry transformation 
W 9 = ±/ 2x2 , there arc Z4 classes of different fermionic 
non-chiral SRE phases in the presence of Z2 x Z 2 sym- 
metry. The nontrivial SPT phases such as [0,1,0] can 
be realized by non-interacting fermions, as shown by the 
coupled wire construction in section |VI| 

Interacting fermionic SPT phases from the 
bosonic SPT phase with Z 2 symmetry 

In the previous discussion of bosonic SPT phases, a 
Z 2 classification was found for bosons with Z2 symme- 
try. Here in a fermionic system with Z 2 x z£ symmetry, 
we can always let the fermions combine to form bosonic 
Cooper pairs which can serve as the fundamental bosons, 
which then form the nontrivial bosonic SPT phase dis- 
Notice that the fermion parity Z 2 
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cussed in section 
symmetry can never be broken and have no effect on the 
Cooper pairs at all. Do these interacting bosonic SPT 
phases lead to an extra Z 2 classification for fermions 
with Z2 x Z2 symmetry, in the presence of deconfincd 
fermions in the low-energy sector? If so these non-trivial 
SPT phases cannot be obtained from perturbing non- 
interacting fermions. However, we show now that the 
bosonic SPT phase with Z2 symmetry is contained within 
the fermion classification discussed previously when there 
are gapless fermions on the edge. And it is a Z2 subgroup 
of the Z4 classes that were found. Thus they can be 
obtained from adding perturbation to a non-interacting 
fermion Hamiltonian. 

Consider one bosonic ^-symmetric SPT state with 
edge variables {<j>\, 2 }, whose symmetry transformations 
are 4> a — > a + 7r, a=l,2 under Z 2 generator g. When 
this state is put together with two fermion Z 2 x z£- 
symmetric SPT states [0,1,0] © [0,1,0] with edge vari- 
ables {(j)R, 0l} and {4>' R , 4>' L }, its edge can be gapped out 
by simultaneously localizing the following bosonic fields 
on the edge: 



'i? -t- <Pl -t- Pi, 

<P'r + 0L + 02, 

1>r- 4>' L - <t>2- 



(107) 



Notice that under Z 2 generator g the edge variables 
{<t>R, <$>l} transform as (j> R -t <p R + w, cf> L -t <p L and 
the same for {<f>' R , <f>'i,}- Notices that the inverse of the 



above bosonic SPT phase is itself, hence we have shown 
that [0, 1, 0] 2 = [0, 1, 0] © [0, 1, 0], i.e. the state obtained 
by putting two fermion SPT phases [0,1,0] together 
is nothing but the bosonic Z 2 -symmetric SPT phase. 
Therefore we conclude that bosonic SPT phase with Z2 
symmetry is contained within the fermionic classification 
Z A . 



2. Gf — Z4, symmetry: Z2 classes 



Also note that when r; = 1 we have the algebra 



9 2 = Pf- 



(108) 



for symmetry group G/ = Z4. Note that Z 2 is a normal 
subgroup of Z4. Since all phases with r\ = 1 are trivial, 
they do not give rise to nontrivial (intrinsic) fermionic 
SPT phases with Z<± symmetry. However as discussed 
before, the bosonic SPT phase of Cooper pairs with 
symmetry Gf/Z( = Z 2 (when fermions are confined) 
always automatically lead to one interacting fermionic 
SPT phase protected by symmetry Gf = Z4. Hence all 
different SRE fermionic phases with symmetry group 
Gf = Z4 have at least a Z 2 classification. These are 
physically related to charge-4e superconductors in two 
dimensions protected by electron charge conservation 
modulo 4. The nontrivial fermionic SPT phase protected 
by Gf = Z4 symmetry cannot be obtained by perturbing 
non-interacting fermions, and therefore they are intrin- 
sic interaction-driven fermion SPT phases (charge-4e 
superconductors in this case) with Z4 symmetry. 

In summary, there are at least Z4 classes of different 
fermionic non-chiral SRE phases in the presence of Z2 x 
Zl symmetry. Bosonic SPT phases with symmetry Z 2 
do not add any new phases. 

On the other hand, in a fermion system with Z4 sym- 
metry (g 2 = Pf or 77 = 1) there are Z 2 classes of different 
fermionic SRE phases. This corresponds to two different 
classes of charge-4e superconductors in 2+1-D. 

Discussion of Results: The fermionic topological 
phases protected by Z 2 x Z^ symmetry form a Z 4 group. 
In comparison, super-cohomology theory 2 -' obtains the 
same number of phases but with group structure Z 2 . An 
advantage of the present formalism is that we can see how 
these phases connect to the bosonic SPT phases with the 
same symmetry, and verify they do not add any new 
phases. 



F. Gf/Zl = Z 2 x Zl symmetry 



The algebraic structure of Z 2 x Z^ x Z 2 group is 



gTgT = e. 



(109) 
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where g is the Z2 generator and time reversal operation 
T is the Z2 generator. In our notation e can be either 
identity element e / for fermions or fermion number parity 
Pf. The associate constraint (27 1 are given by 



(W 9 ) 2 = {W T f = (W 9 W T ) 2 = J 2x2 , 
(W 9 ) T KW 9 = ~(W T ) T KW T = K, 



(110) 



(1 + W 9 )8<j) 9 = r/ g n 



(1 - W T )8<}> T = r, T TT 



(1- W T W 9 ){5cf> T - W 1 



') = r/TT 



(111) 



where 77, t/t, i] g = 0, 1. W e can always choose a gauge so 
that W T = (j £ and from (fllbl IF 9 = ±J 2x2 . We haven't 



found a microscopic realization of symmetry transforma- 
tion W 9 = —^2x2 so far, therefore we put the discussions 
of W 9 = —^2x2 case to the Appendix. Here we'll focus 
on symmetry transformation W|x2 case. 



For W 9 = I2 x 2 the inequivalent solutions to ( 111 I are 




t* T = K j , a,r],r) T ,r] g = 0, 1. 



If r] g + ri — tit = the variable </>i + 02 can be localized 
without breaking any symmetry. If 77 = 77X = the 
variable <pi — 2 can be localized without breaking any 
symmetry. 

Note that when rj g = we have g 2 = e/ and hence the 



symmetry group is G/ = Z2 x Z2 x Z2 with T 
if 77T = or T 2 = P f if ? ?T = 1. When ?7 g = 1 



: e f 
the 



other hand we have g = Pf, and symmetry group for 
fermions is Gf = Z 4 x or G/ = Z4 x Zj. 



1. G / 



■ Z2 x x Zjj symmetry: Z4 x Z2 classes 



There are 4 nontrivial SPT phases with r/ g = 0: 
they have r\ — t\t — 1 mod 2 and a = 0,1. Let's 



label a state with symmetry transformations (1121 as 
[Vg 7 Vt >i]>oi\' We already showed [1,0,0, a] — [1 , 77 + 
1,7/, a] = [0,77, 77, a] = e Z2xZ T xZ /. When a [0,77 + 1,77,0] 
state with edge variables {<pi, (f>2 } is put together with a 
[0,77' + 1,77', a]" 1 state with edge variables {(/)[, <fi' 2 }, its 
edge cannot be gapped without breaking the symmetry 
by localizing independent bosonic fields </>2+(/> 2 }- 
Therefore we know [0,1,0, a] = [0,0,1, a] is the same 
nontrivial SPT phase. In the case ?7 g — = t)t and 
77 = 1 the algebra of Z 2 x Zj x Z% group is 



As discussed in Ref. HT] this is the same as rj g = = 77 
and 77T = lj since one can always redefine the anti- 
unitary time reversal as T' = gT. Just as discussed 
earlier for fermionic SPT phases with Z2 x symmetry, 
it's easy to verify that [0, 1, 0, 0] = [0, 1, 0, 1] _1 and when 
four [0, 1,0,0] states are put together their edges can be 
gapped without breaking any symmetry i.e. [0, 1, 0, 0] = 
e z 2 xz T xz f ■ ^ S a resu lt [0,1,0,0]™, n = 1,2,3 are the 
only three nontrivial SPT phase, whose edge cannot be 
gapped without breaking the symmetry. Hence all dif- 
ferent phases with W 9 = I 2 x2 form a Z4 group for 
Z2 x Z2 x Z 2 symmetry (?7 g = 0). 

For the same reason mentioned earlier for fermions 
with Z2 x Z^ symmetry, here in the presence of Z2 x x 
Z 2 symmetry, we can obtain interacting SPT phases from 
bosonic SPT phases of fermion pairs with Z2 x sym- 
metry. Note there are classes of bosonic non-chiral 
SRE phases with with Z2 x Z^ symmetry. Do these lead 
to an extra Z| group structure for fermions with this 
symmetry? As before, we now show that this is not the 
case. These phases are already accounted for within the 
fermionic classification when fermions are present on the 
edge. While one of these Z2 classes is contained within 
the Z4 classification, the other is connected to the triv- 
ial class, in the presence of fermions. Hence it brings 
in an extra Z2 structure due to the bosonic SPT phase 
of strongly-bound Cooper pairs where fermions are ab- 
sent (confined) on the edge. Such a Z2 index cannot be 
obtained from perturbing free fermions. 

Again let us consider a bosonic Z2 x Zj-symmetric 
SPT state [l,7i2, n]b with edge variables {tjji,^}, which 
transforms as 




under Z2 generator g and time reversal T . When this 
state is put together with two fermion Z2 x Z£ -symmetric 
SPT states [0,1,0,0]/ ® [0,1,0,0]/ with edge variables 
{(^r^l} and {4>'r,4>'l}' edge can be ga pped out by 
simultaneously localizing the bosonic fields (107) on the 
edge if n = 1. This means 

[1, n 2 , l] b © [0, 1, 0, 0]/ © [0, 1, 0, 0]/ = e Z2XZ r xZ f . 

Similarly one can show that 



[1, 77 2 , o] h e [0,1,0,0]/® [0,1,0,1]/ = 

= [0,1, 0,0]/ ©[0,1, 0,1]/. 



J Z 2 xZ,' 



gTgT = Pj 



(113) 



again by localizing the same bosonic variables (1071. 
Hence we have shown that [l,7i2,l]b = [0,1,0,0]/® 
[0,1,0,0]/ and [l,n 2 ,0] 6 = [0, 1, 0, 0]/ © [0, 1, 0, 1]/. So 
in the presence of fermions, all the bosonic SPT phases 
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protected by Z 2 x Zj x z| can be obtained by perturb- 
ing non-interacting fermions (constructed by putting sev- 
eral non-interacting [0, 1, 0, 0]/ together and add interac- 
tions). Therefore we conclude that bosonic SPT phase 
with Z 2 x Zj symmetry can all be obtained from perturb- 
ing non-interacting fermions in the presence of Z2 x Zj 
symmetry. However when fermions are confined in the 
low-energy sector (on the edge), the bosonic SPT phase 
[1,1, 0]b of strongly-bound Cooper pairs indeed bring in 
new fermionic SPT phase for fermions. Hence in total 
there are Z 4 xZ 2 classes of different fermionic (non-chiral) 
SRE phases with Z2 x Zj x Z 2 symmetry: Z 4 can be ob- 
tained from perturbing free fermions, and the other Z 2 
is from bosonic [1, ri2, 0]b phases of molecule-like Cooper 
pairs protected by Z2 x Zj symmetry where fermions are 
confined in the low-energy limit. The last Z 2 cannot be 
obtained from perturbing free fermions. 

Discussion of Results: While this particular sym- 
metry class cannot be discussed within group super- 
cohomology theory^, recent wortP^^U have approached 
this problem from another angle, by starting with non 
interacting fermions (which have a Z classification with 
this symmetry) and then turning on interactions. They 
find a Z§ classification that survives interactions. Odd in- 
teger members of this series have an odd number of pair 
of Majorana modes at the edge, that move in opposite 
directions. Although apparently quire different, these re- 
sults are consistent with ours due to the following. Since 
we are unable to deal with unpaired Majorana modes, 
only the even members of the series are captured (hence 
Z 4 in our Z 4 x Z 2 classes here). An advantage though 
is that this classification of topological phases which are 
stable to interaction emerges directly from the formal- 
ism, without the need to begin from free fermions. And 
our formulation also includes a Z 2 classes that cannot be 
obtained by perturbing free fermions. 



2. Gf — Z4 x Z2 symmetry: Z| classes 

On the other hand, when r\ g — 1 we have g = Pf = 
(—l) f and hence the corresponding symmetry group is 
actually Z4 x Zj, where Z 2 is a subgroup of Z 4 . In 
this case the nontrivial SPT phases have rj g — 1 and 
V = Vt = 1, a = 0, 1. The algebraic structure of the 
symmetry group Gf is given by 



g 2 = T 2 =gTgT = P f . 



(114) 



and hence 
2 



the 



It's easy to check that Tg = g 

symmetry group is actually Z 4 x Z 2 with T z = Pf. 
We still label the phases with symmetry transforma- 
tions (1121 as [i] g , t]t, i], a}. When a [1, 1, 1,0] state with 
edge variables {</>i, </> 2 } is put together with a [1, 1, 1, 1] _1 
state with edge variables {0'i,0' 2 }, its edge cannot be 
gapped without breaking the symmetry by localizing in- 



dependent bosonic fields {cf>i 



,cf>2 + </> 2 }. Therefore 



Now let's put two [1,1,1,0] states together with edge 
variables {4>l, 4>r} an d {4>'l> 0i?Jj it's easy to see the edge 
states will be all gapped out by simultaneously localiz- 
ing the following bosonic variables: {<fi 7? — <j)' L , <j)' R — 4>l}- 
Hence we have [1,1, 1,0] 2 = e Z4xZ T and these differ- 
ent intrinsic fermionic SRE phases with Z 4 x Z-f form 
a Z 2 group. The nontrivial fermionic SPT phase can 
be obtained from free fermion band structures, just like 
Gf = Z 2 x Zj x Z{ case. 

In the case of a fermion system with Z 4 x Zj sym- 
metry, the corresponding bosonic system of Cooper pairs 
have Z2 x Z-f symmetry and has a Z 2 topological classi- 
fication. For fermions in the presence of Z 4 x Zj symme- 
try, these bosonic SPT phases of strongly-bound Cooper 
pairs (with Z 2 x Zj symmetry and a Z 2 classification) 
will bring in new phases, in addition to the the non- 
trivial fermionic SPT phase [1,1,1,0] (with a Z 2 group 
structure). In these new phases the fermions can be ei- 
ther confined in the low-energy sector (associated with 
[1, 1, 1, 0] 2 state) or deconfined (associated with [1, 1, 1, 0] 
state). Hence all different fermionic SRE phases with 
Z4 x Zj symmetry have a Z 2 classification, where one 
Z 2 structure comes from bosonic SPT phases [l,n 2 ,?x]b 
of Cooper pairs, and the other Z 2 associated with fermion 
state [1,1,1,0] are intrinsic properties of fermionic sys- 
tems. 



Gf = Z4 xi Z% symmetry: 



Z| classes 



In this case we have r\ g = \ and r\ — and therefore 
9 



I'. gTgT = e f , (115) 
T 2 = e f or P f . 



In this case there are no intrinsic fermionic SPT phases 



with G 



f 



Z2 symmetry. 



However bosonis SPT 



phases of Cooper pairs with G f jZ^ — Z 2 x Zj symme- 
try leads to Z 2 classes of different fermion SRE phases. 
Hence there are at least Z 2 classes of different fermion 
non-chiral SRE phases with Gf = Z 4 x Zj symme- 
try. These Z 2 classes of phases physically correspond 
to different charge-4e superconductors with time rever- 



sal symmetry T 



Recently, the 



[1,1,1,0] = [1,1,1,1] is the same nontrivial SPT phase. 



possibility of realizing charge-4e superconductivity (four 
fermion condensates), in imbalanced cold atomic gases^ 
(which break time reversal symmetry) and also in cer- 
tain cuprate superconductors 50 (which preserve time re- 
versal) has been discussed. While these phases were non- 
topological, the prospects for realizing topological phases 
with these symmetries remains to be seen. 



VI. COUPLED WIRE CONSTRUCTION OF 
BOSONIC AND FERMIONIC SPT PHASES 

In the previous sections we showed how to classify dif- 
ferent SPT phases in the K matrix + Higgs formulation. 
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spin 

charge 




Figure 2: Schematic illustration of interwire coupling terms 
which stabilize the bosonic SPT phases protected by U(l) 
symmetry, with Hall conductance a xy = 2q. Solid horizon- 
tal lines stand for quantum wires of charged bosons (each 
carries unit U(l) charge) while dashed horizontal lines repre- 
sent quantum wires composed of neutral (say spin) degrees of 
freedom. Dashed and solid arrows illustrate the two interwire 
coupling terms in ( |130[ ) that gap the bulk, but leave behind 
nontrivial edge states. 



The edge structure of SPT phases is explicit in this for- 
mulation: e.g. the edge of a bosonic SPT phase is char- 



acterized by bare action ( 34 ) , Kac-Moody algebra ( 35 1 



as well as symmetry transformation rules ( 22 ) for uni- 
tary symmetry g and ( 24 1 for anti-unitary symmetry h 
on the bosonic variable {0 a , a = 1, 2}. However, a more 
microscopic construction of these 2 + 1-D SPT phases 
is still lacking in this formulation. In this section we 
present a microscopic construction of these SPT phases 
in the anisotropic (quasi-lD) limit, from an array of cou- 
pled one-dimensional quantum wires. This approach has 
been used to construct AbeliarP^ and non-Abelian FQH 
states^. We first give a short introduction to the cou- 
pled wire construction, and then use the this method to 
explicitly construct bosonic SPT phases in the presence 
of symmetry group G — U(l) and G = U(l) x Z%, as well 
as fermionic SPT phases with symmetry group Zi x Z% '. 
Generalizations to other symmetry groups are straight- 
forward. 



A. Coupled wire construction in a nutshell 

Consider an array of uncoupled identical one- 
dimensional quantum wires, each wire being described 
by a non-chiral Luttinger liquid. The bosonic fields 
associated with the Luttinger liquid in the l-th wire 
(1 < I < N w , N w being the total number of quantum 
wires) are {6i(x), tpi(x)} satisfying the following commu- 
tation relation 



[O m {x),<pi(y)] 



i — Sign(x - 



(116) 



<p{x) is a bosonic phase field, while 9(x) describes the 
density fluctuations in the Luttinger liquid. The long- 
wavelength density fluctuations on l-th wire is given by 



Pi(x) - pi= d x 0i(x)/ir. 



(117) 



where pi is the average particle (boson or fermion) den- 
sity in the l-th wire. In terms of these two variables, 
the Luttinger liquid Hamiltonian for decoupled quantum 
wires is given by 



2tt 



da; 



{d x 0if + gi{d<pif . (118) 



where gi > is the Luttinger parameter for the l-th wire. 

The idea of coupled wire construction is to add inter- 
wire and intra-wire interactions between electrons as well 
as tunneling between wires. For example the forward 
scattering terms between different wires is written as 



Hfc = Y,I dx (d<p k ,d6 k )M 




(119) 



where Mj^ are all 2 x 2 matrices describing the forward 
scattering interactions between the j-th and fc-th wire. 
Other inter-channel scattering terms in general have the 
following form: 



cos 



i J2i (miOi(x) + nupi(x)) + OL{ muni } 



(120) 



where a{ mijT!! } are real constants and {mi, n{\ all take in- 
teger values. For a bosonic system in the absence of any 
symmetry (or associated conservation laws), the above 
interwire scattering term must satisfy the following con- 
dition: 



mi = mod 2, V 1 < I < N v 



(121) 



This is because the boson density fluctuations are mainly 
contributed by density waves at vector q n ~ 2itp~in, n € 
Z and the density fluctuations at q n is given by 



o l n (x) cx e i»(a*ft*+iWi(*)) i ne z, l<l< 



N„ 
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for the Z-th quantum wire. For a fermionic system on 
the other hand, in the absence of any symmetry the con- 



straint on interwire scattering term (120) is 



mi = ni 



mod 2, V 1 < I < N„, 



(122) 



This is because inter channel scattering terms must be 
composed of single fermion operators: i.e. left mover 
ipl ~ exp[i((/3; — 9i — irpix)] and right mover ipf' ~ 
cxp[i(<y3; + 6i + wpix)]. The presence of symmetry group 
G will lead to further constraints on inter-wire cou- 



pling terms (120) and forward scattering terms (119 1 



symmetry-allowed scattering terms must transform triv- 
ially under any symmetry operation. The bare Luttinger 
liquid Hamiltonian together with symmetry-allowed for- 
ward scattering (1191 and inter channel scattering (120) 



forms the generic Hamiltonian for a coupled wire con- 
struction: 



H = H 



LL 



H 



FC 



(123) 



In the coupled wire construction (123 1, one can prop- 



erly choose a set of symmet ry-al lowed interwire scatter- 
ing terms {0^- m '' n ^{x)} in (fijl, so that 



^ mini = ^(m/nj + mjri;) = ^ 



m{n{=0. (124) 



for any two terms { - mi ^(x) and 0^ m 'i> n '^{x) in Hamil- 
tonian (1231. Therefore the set of bosonic variables 
{^2i(rni0i + nitpi)} can be simultaneously localized at 
classical values by minimizing the interwire scattering 
terms {0^ m '' n '^ (x)}. When chosen properly all degree 
of freedom in the bulk will be gapped by these interwire 
scattering terms, and the only low-energy degrees of free- 
dom left free are on the edge. To be specific the variable 
V\0\ + ?i</?i on the left edge I = 1 will remain gapless if 



•p\Ti\ + qirrii — 0, V {mi, ni} in ( 123 ). (125) 



Note that one can always tune the forward scat- 
tering terms ( |119 1 so that interwire coupling terms 
are relevant in the renormalization group 
sense. Then one expects the coupled wire system will 
be driven into strong coupling phase of interwire scat- 
tering {0^" ll ' ni ^ (x)}, and hence all bosonic variables 
{^2i{mi0i + nitpi)} will be localized simultaneously at 
classical values. 

One of the simplest example is the Laughlin state 77 
of spinless fermions, i.e. FQH state at filling fraction 
v = 1/m, m =odd integer. The interwire scattering 
terms whose strong coupling phase corresponding to the 
Laughlin state are^ 



I m 



N w -1 

E 

i=i 



dxCi cos 



<Pi - w+i - m(9i + i+ i) 



Its gapless variable on the left edge is 4>^ m = <fii 4 
which satisfies the following Kac-Moody algebra^ 

[ d x<t>i/ m (%), 9 v <Pi /m (y)] = 2irmid x S(x - y). 



It's easy to verify that </>f/ m satisfies condition (125| and 
is the only gapless degree of freedom on the edge 



B. Bosonic SPT phases with U(l) symmetry 



As discussed in section IV in the presence of U(l) 



symmetry there are Z (integer group) different classes of 
bosonic non-chiral SRE states, which are labeled by their 
U(l) charge vector t = (l,q) T . The bosonic variables 



{^1,^2} on its edge satisfies Kac-Moody algebra (351. 
Under group element of symmetry group U{1) they 
transform as 




(126) 



The nontrivial SPT phases correspond to different in- 
tegers q 7^ 0, whose edge cannot be gapped out without 
breaking the U(l) symmetry. Here we present an explicit 
construction of these SPT phases with U(l) symmetry in 
the coupled wire approach. 

We start from an array of quantum wires (1 < I < N w ) 
where each wire is composed of two chains: a chain of 
charged bosons (each boson carries a unit of U(l) electric 
charge) and a spin chain. Each chain forms a c — 1 Lut- 
tinger liquid described by bosonic fields: {tpi(x), 9[(x)} 
for the spin chain, {ipi (x) , df (x)} for the chain of charged 
bosons in the Z-th wire. These bosonic fields satisfies the 
commutation relation (|116|): 



K(*),<P?(y)] =i^Sign(x 



(127) 



where a, j3 = c/s denotes charge/spin degree of freedom 
and 1 < m, I < N w are the wire index, ip^ are phase fields 
of charged bosons while 2d x 8f describes charged boson 
density fluctuations. For the spin chain d x 0[(x) ~ Sf(x) 
and exp[iip{(x)] ~ S^(x). Without interwire scattering 
terms, the bare Hamiltonian density of the system takes 
the form (|118|) of Luttinger liquids 



n 



LL 



EE 

c/ 5 



1 = 1 



2tt 



9\ 



-(Mf) a +flf(ty>P) a • (128) 



The ?7(1) symmetry associated with {ipf(x),6f(x)} bo- 
son charge conservation leads to the following symmetry 
transformations for the bosonic fields: 

V, c (aO -► UA*<tf{x)UZl = <pf(x) + A0, (129) 
U A<i> = e iA 'l'f dx 'Ei 2d * s '( x \ O<A0<2tt. 

The other fields Of, ipf, Of are invariant under the above 
U(l) charge rotation U A <p. 

In the presence of the above U(l) symmetry associated 
with boson charge conservation, the different phases la- 
beled by charge vector t = (l,q) T are stablized by the 
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following inter- wire coupling terms: 



«(i,«ri = S^i" 1 Q cos(^ - - 2Q\ + *0 + Di 



L (l,-3) 



cos (ipf - <pf +l + q(<pf - <p c l+1 ) - 2(6f +1 



'1+ 



i) + A0 



(130) 



{^x(af), 0i (a;)} are gapped for the same reason while vari- 
ables {(pj(x), <pf(x)} are left gapless. For the I = N w wire 
on the right edge, things happen in the opposite way: 
variables {fa N (x),^ (x)} are gapped while variables 



{<t>l 



(x)} remain gapless. 



where C;,-D/,A;,AJ are real constants. A pictorial illus- 
tration of the above interwire scattering terms is given 
in FIG. [2] Clearly the above interwire scattering terms 
all satisfy constraint ( |121[ ) for bosonic s ystem s, and they 
are also invariant under U(l) rotation (1291. 



As argued in Ref. 151)521 one can always choose proper 
forward scattering terms (119 1 to make the above inter- 



wire coupling terms become relevant and drive the system 
into their strong coupling phase. Notice that the argu- 
ments of the above cos terms commute with each other, 
so they can be localized at certain classical values simul- 
taneously. It's straightforward to show that all bosonic 
fields in the bulk with 2 < I < N w — 1 are gapped while 
the gapless edge states on the left edge / = 1 are described 
by variables {<f)\(x), <pf{x)} defined as 



# = ft 



# = f°v 



(131) 



They transform exactly like {fa, fa} in (1261 under 



charge U(l) symmetry (129 1. Besides they also obtain 
the Kac-Moody algebra ( 35 ) for a bosonic non-chiral SRE 



system. As a result the strong coupling phase of interwire 
couplings (130) is nothing but the bosonic SPT phases 
labeled by charge vector t = (l,q) T with charge U(l) 
symmetry. 

Now let us elaborate on why the interwire coupling 
( 130[ ) can gap out everything in the bulk and leave 
variables (131) on the edge. In addition to variables 
{(pi (x), 4>i(x)\ in (1311 one can define another pairs of 
variables {<pj (x) , <pf (x)} as 



fa = <pf- 26, 



1 ' 



■<m- 



(132) 



They also satisfy Kac-Moody algebra ( 35 1 except for an 



extra minus sign for all commutators. Notice that the two 
pairs of variables {fa (x) , fa (x)} and {fa(x),fa(x)} com- 
mute with each other. They are just a linear combination 
of the original charge and spin variables {ipf, Of, ipf , 6? }. 
The interwire scattering terms ( 130 1 can be written as 



(1,5) 



— 1^,1 = 1 



Ci cos( 



+£>, cos (fa - cf>t +1 + a;) 



C. Bosonic SPT phase with (7(1) XI Z% symmetry 



As discussed in section IV in the presence of U(l) xi 
symmetry there are Z2 different classes of non-chiral 
bosonic SRE phases in 2 + 1-D. Among them there is only 
one nontrivial SPT phase, whose edge cannot be gapped 
without breaking the U(l) x Zj symmetry. Its gapless 
edge is described by bosonic fields {(j>\, fa} which satisfies 
the Kac-Moody algebra (35 1. Under charge U (1) rotation 
the two bosonic variables transform as ( |126 1 while under 
time reversal T they transform as 



T : 




(133) 



for the nontrivial SPT phase. 

Such a SPT phase is nothing but the strong coupling 
phase of interwire scattering terms (1301 with q = 0. Its 



gapless bosonic fields on the left edge / = 1 are 



(134) 



<t>\ = <p\ + 26? 



while on the right edge / 
are 



N w the gapless boson fields 



fN„ 
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N„ 



(135) 



fN m - 

Under time reversal T the original boson fields 
{cp?, 9f, iff, Of} naturally transform as 



= ft +7T, 



T0fT~ 
T6fT~ 



'I • 



since all components of the spin S z ~ d x 9 s , S + ~ 
exp[i(/3 s ] change sign under time reversal. It's easy to 
verify that interwire scattering terms ( 130 1 with q = 



is invariant under time reversal as long as we choose 
A; =0. Hence these interwire couplings are allowed 
by symmetry. It's also straightforward to show that the 
pair of bosonic fields, i.e. both {<fij (x) , <f>f (x) , I = 1} and 
l<t>i(x ),(t>i( x), I — N w } transform in the same way as 
(1261 and (1331 under U(l) X symmetry. Hence the 
strong coupling phase of interwire couplings (1301 with 



q = indeed corresponds to the nontrivial bosonic SPT 
phase in the presence of U(l) X Zl^ symmetry. 



Hence everything in the bulk i.e. {fa' (x),fa' (x), 2 < 
I < N w — 1} are all gapped since they don't com- 
mute with at least one of interwire scattering terms in 
(1301. For the I — 1 wire on the left edge, variables 



D. Fermionic SPT phases with Z2 X Z% symmetry 

Here we show that fermionic SPT phase [77 = 0, t\ = 
1, t% = 0] with W 9 = ^2x2 discussed in section [V] can be 
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constructed in the coupled wire approach. Its coupled 
wire construction also indicates this SPT phase responsi- 
ble for the Z 4 group structure can be obtained from non- 
interacting fermion band structures. The edge variables 
{<fii, ^2 } satisfies Kac-Moody algebra (73 1 and transform 
in the following way under Z2 symmetry g: 



(136) 




and under fermion parity Z 2 = {e, Pf}: 




(137) 



Consider right now each quantum wire contains elec- 
trons of both spins i. e. two left movers V'f/^ an d two right 
movers ipy^'- 



exp 
' exp 



+ k 



t/i x ) 
k r/i x ) 



where bosonic fields (pi_ a and 0/ !<T satisfy commutation 
Let's assume here Z2 symmetry g is the 
iber parity of spin-f fermions, which is nat- 



relation (116 



fermion num: 
urally realized by 

g : ipi,f(x) (fi t f(x) + ir. (138) 

while 0i^,Lpi ^,9i^ remain invariant under g. On the 
other hand under total fermion parity Pf we have 

P f : <Pl,1-/i,i x ) -> <PlA/l( x ) + n - ( 139 ) 

where 9i^,6i t i remain invariant under Pf. 
By defining the following variables 



i ' — Via + "i,ti 



$ = VIA 



and 



*f = VIA 



t>i=m + i,i- 



(140) 



(141) 



It's easy to verify that both {^>; , 4>f} and {(frf 1 , 0f } sat- 
isfy Kac-Moody algebra ( 73 I and the symmetry trans- 
formations (136) and (137). The two pairs of variables 



commute with each other. The two variables {0f-,</>; L } 



are nothing but the right mover for spin-f and left mover 
for spin- J,. Clearly the following single-fermion tunneling 
terms 



^z 2 xz| 



+ + A cos( 



AO 



will gap out everything in the bulk in its strong coupling 
phase. The gapless variables a re \4 >x, <Pi } on the left 
edge 1 = 1: they do transform as ( 136 1 and ( 137 ) . On the 



VII. CONCLUDING REMARKS 

We have discussed an algebraic method to systemati- 
cally classify interacting topological phases in two dimen- 
sions in the absence of topological order. The key devel- 
opment is a general formalism for incorporating symme- 
try transformations into the K matrix formalism. Vari- 
ous examples of interacting boson and fermion topologi- 
cal phases that are well defined in the presence of disorder 
were presented. The method provides both long wave- 
length information of these phases (bulk effective field 
theory and edge theory) as well as suggests microscopic 
realizations in model systems (such as quasi ID realiza- 
tions presented here). It also opens the door to study 
various transitions out of these phases - e.g. topology or 
symmetry changing transitions driven by disorder and 
interactions. Future work will focus on extending these 
results to symmetry protected distinctions between topo- 
logically ordered states, and extending this formalism to 
d = 1 and d = 3. It remains to be seen if a more gen- 
eral formalism that subsumes the present one can be de- 
vised, which, for example, can handle unpaired Majorana 
edge modes. A deeper understanding of the remarkable 
connection between this formalism and the Borel group 
cohomology/super-cohomology is also needed. Perhaps 
the most important question is to determine how topo- 
logical phases of the interacting variety can be obtained 
in an experimentally realistic setting. We leave these 
questions for future work. 

While completing this manuscript Ref. 59 appeared 
which studies the specific case of G = U(l) x Z 2 (topolog- 
ical insulators) using a K matrix approach. Our results 
agree in the areas where they overlap. 
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Appendix A: A short note on GL(N, Z) 

GL(N, Z) is the group of all N x N unimodular ma- 
trices. All GL(N, Z) matrices can be generated by the 
following basic transformations (i =/= j): 



right edge I = N w the gapless variables are ,<p^ }. 
Since the interwire scattering term which stablizes this 
SPT phase is just a single-electron tunneling term, we 
expect such a fermionic SPT phase should be realized in 
a non-interacting band structure with symmetry Z2 x Z 2 . 



Sab = Sa,b(l ~ S a< i)(l ~ S a< j) + 5 a jSb,i ~ Saa<>bj: 



D a .b — S a ,b — 2<Ljv#, 



b.N- 



(Al) 



TW'K will add the j-th row of matrix K to the i-th 
row of K, while S^'^K will exchange the i-th and j-th 
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row of K with a factor of —1 multiplied on the i-th row. 
£>K will just multiply the iV-th row of K by a factor 
of -1. KTf'J), K5' tJ > and KD correspond to similar 
operations to columns (instead of rows). A subgroup of 
GL(N, Z) with determinant +1 is called SL(N,Z) and 
it's generated by {T^\ S^}. 

As a simple example when N = 2, group GL(2,Z) is 
generated by the following basic transformations: 



S = 



-1 



T = 



1 1 



v 1 / V° v 

The following results will be useful 



1 n 
1, 



(STS) n 




Appendix B: A theorem on 2 x 2 K matrices with 

det K = -n 2 

In general a gauge transformation which relabels the 
quasiparticles in the K matrix formulation ([2| is imple- 
mented by a GL(N,Z) matrix W as shown in ^ and 
^ . Therefore in the absence of any symmetry, any two 
K matrices related by (|8| are equivalent to each other 
i.e. 



K ~ M T KM, V M e GL(n, Z). 



(Bl) 



We use ~ to denote the equivalency. In the presence of 
U(l) symmetry with a charge vector t the equivalency 
requires 

{K,t} ~ {M T KM, M T t}, V M e GL(n, Z). (B2) 

In the absence of any symmetry, here we prove the 
following theorem: any 2 x 2 K matrix with determinant 
det K = — n 2 can be transformed into the standard form 



n 



, 0< a< 2n- 1. 



In the special case of SRE phases with n = 1 
fo n\ 



\n 2a I 

'o n 
n 2a + 1 



< a < n — 1 for bosons, 
, < a < n — 1 for fermions. (B3) 



A generic 2 x 2 K matrix with determinant — n 2 can 
be written as 



K 



2x2 



a n + k \ 
n + k b I 

ab = k(2n + k), a,b,k € Z. 



(B4) 



Apparently for a bosonic system a, b are both even in- 
tegers and fc is also an even integer. For a fermionic 
system there are two possibilities: a, b, k are all odd inte- 
gers; k — or — 2n and one of a, b equals zero while the 
other being an odd integer. 

Notice that under GL(2, Z) transformations a x and 
ia y (cr a , a = x,y, z are Pauli matrices) we have 



(B5) 




If a = or b = we have k = or k 



-2n i.e. n + k 



±n in (B4|. Using relations (B5 1 and generator T in ( A2 1 



one can show 



K 



n 




and (B3l can be easily verified. 



If ab 7^ 0, without loss of generality we can assume 
that | a | < \b\ and therefore \a\ < max(|fc|, \2n + k\) — 
max(|Ki 2— n\, |Ki 5 2+^|)- We use the following strategy: 

1 — sign(a(2n + k)) 







1 



if |fc| < \2n + k\ choose M = 

+ n| |Ki, 2 + nl-lal; if \k\ > \2n + k\ 
sign(afc) 



(Bl) so that 



choose M = 



Kl,2 - 

1 





1 



in (Bl) so that |Ki, 2 - 



|Ki ( 2 — n\ — |o|. The value of max(|Ki.2 — n\, \K.i,2 + 
n|) will decrease monotonically when this procedure is 
repeated and finally one will end up with a 2 x 2 K matrix 
whose off-diagonal elements are ±n. This means ab = 
in (B4|. Therefore theorem (B3l is proved. 



Appendix C: A theorem on bosonic 2 x 2 K matrices 
with detK = — 1 and charge vector t 

In this section we prove the following theorem (we use 
(a, b) to denote the greatest common divisor of two inte- 
gers a and 6): for a 2 + 1-D bosonic system any K with 

detK = — 1 and a charge vector t = I 1 ] with (ii,£2) = 

V 2 / 

1 is equivalent to K = | ^ M and t=| ^ | , I £ Z k 
a GL(2,Z) gauge transformation, i.e. 



{K,t 




if(ti,t 2 ) = l, 



First we notice that according to the Euclidean division 
algorithm on integers Z, for any pairs of integers e.g. t\ 
and ti here, there is always a list [qi, g 2 , ■ • ■ , g n +i] and 
[ri,r 2 , • • • ,r n ] such that (let's assume |fi| > \t2\ without 
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loss of generality) 



by the following GL(4, Z) transformations: 



h = Qih + ri, 
h = q2n + r 2 , 
ri = q^2 + r 3 , 

r n-2 = 1n T n-\ + r m 

r n -i = q n +ir n - 

where r n — (ii,i 2 ) is the greatest common divisor of t\ 
and £2; and 1 < |r m +i| < \r m \, Vto. Therefore one can 
always find two integers u± and 112 such that 



(*i,t 2 ) = r„ = r„_ 2 - ft^n-l = • 
As a result for (ii,i 2 ) = 1 we have 



*1«2 - *2«1-(C2) 



and hence 



T _ 



h Ul e GL(2,ZX.C3) 

,*2 "2/ 



{K,t} ~ {K 7 = (Mo^yKMo -1 , f * 



(C4) 



as long as (ti,ta) = 1- In the following we prove that an 
arbitrary 2 x 2 K matrix with det K = — 1 for a bosonic 

system is equivalent to the standard form ^ by a 

gauge transformation which keeps the charge vector 

invariant. To prove this we need to enlarge the Hilbert 
space by introducing a 4 x 4 matrix K = K2X2 © o~ x and 
associated charge vector t = (1,0, 0, 0) T . This describes 
a direct product of the original 2x2 bosonic K matrix 
with a U(l) charge conservation and another trivial 2 x 
2 bosonic K matrix without any symmetry. A generic 

form for K is [ ^° 2fc + 1 \ ga |.j g £„ m „ a b — k(k + I) 
\2k + l 2b J y ' 

(detK = — 1). One can prove that 



{K 



/ 


2a 


2fc + l 








2k + 


1 2b 
























1 


1 





\ 








1 






w 




(0 


1 0) 




(A 








1 


2/0 







}. 

















1 













1 oj 




W 







(C5) 



/0 1 
1 2b 

2{ab-k) 
\0 -26 

fo 1 o\ 

1 26 a' 6' 
a' 1 
\0 6' 1 0/ 



/2a 2k + 1 0\ 
2fc + 1 2b 0,. 

IVLl 

1 
V 10/ 



2(ab-k) -2b 

2a(ab-2k) 2(k - ab) 

2(k -ab) + l 2b 

[0 1 



(M 2 T 



1 2(b-a'V) 

\0 

a, 6, fc, a', b' € Z. 



1 

/ 

o\ 



1 

1 0/ 



M, 



where we defined 



Mi 



M 2 




1 


2k 

(l 



- 

\o - 




a 
1 

-2ka 2k 

o\ 

1 
-b' 1 
-a' I J 



\ 

-1 





V 



dot 



- 1 . Notice that the 



We have used the theorem ( B3 1 for a n = 1 
bosonic system since one can easily verify that 

2a(a6-2£;) 2(k-ab) + l 
2(k -ab) + l 2b 
charge vector t = (1,0,0,0) T remains invariant under 
these GL(4,Z) transformations. Combing relations (C4| 
and (|C5| we have proved theorem (CI I. 



Appendix D: Faithful vs. unfaithful representations 
of the symmetry group 

Transformations {W 9a , <5</>f a } form a faithful or un- 
faithful representation of symmetry group G: By solv- 
ing the constraint equations ( 26 I and choosing a proper 
gauge in ( p9| , one can obtain a set of transformation 
rules {W 9a , 54> 9 j a } as the solution. Apparently the trans- 
formations {W 9a jScf) 9 . 11 } always form a representation of 
the symmetry group G of the system (or the Hamiltonian 
of the system) in the sense that 

y 91,92 eG: W 9192 = rj^W^W 91 , (Dl) 

(5051-32 = £032 +m W 92 S(j) 91 . 

where 771 = ±1 if g\ is a unitary (anti-unitary) symme- 
try and r\i is associated with g 2 - This representation of 
group G is faithful if and only if the identity element 
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e is the only symmetry element under which all bosonic 
quasiparticle fields {^Zj '/</>/} on the edge (or {H/^/ 1 } 
in the bulk) transform trivially. In other words, under 
any element g ^ e of symmetry group G at least one 
bosonic quasiparticles l T <f) satisfying (111 will transform 
nontrivially. 

In contrast to faithful representations, an unfaithful 
representation {W 9a , 54> 9 j a } of symmetry group G means 
there exists a nontrivial subgroup G^ of G, so that under 
any symmetry element g G G^, all bosonic quasiparticle 
fields li<Pi} on the edge (or {Jl/^/} * n the bulk) 
transform trivially. This means for a phase described by 
K matrix and symmetry transformations {W 9a ,<5</>j°}, 
its edge states can be gapped by condensing the bosonic 
quasiparticles without breaking the subgroup G^, of sym- 
metry group G since under any symmetry g £ G^ all 
bosonic quasiparticles are left invariant. As a result when 
the edge is gapped, the symmetry group G of the Hamil- 
tonian breaks down to its subgroup, the ground state 
symmetry group G^. As a result the symmetry break- 
ing phases can be naturally incorporated in the K ma- 
trix+Higgs formulation. 



Apparently there is only one (Zi class) trivial phase 
e u(i)xz T with U(l) x Z2 symmetry, whose edge states 
can be gapped without breaking the symmetry. 



2. Zn x Z 2 symmetry 



The algebraic structure for Zn x group is given by 



g N = T 2 = Tg^Tg = e. 



(E6) 



where g is the Zjy symmetry generator and T is time re- 
versal. The associated constraint equations for symmetry 
transformations are 



W g W T {W g )- x W T = J 2x2 , 
(7 2x2 + W 9 W T (W 9 ) 



(E7) 
(E8) 



W 9 (l - W T (W a )- 1 )S(f) T = I J mod 2tt. 



in addition to p7p8| and (|53p4 1 



Appendix E: Other Bosonic SPT Phases: 

1. U(l) x Z2 symmetry: Zi class 

In contrast to the U(l) x symmetry discussed in the 
previous section, here we study a direct product of U(l) 
and time reversal Z^ symmetry. This can be realized 
by time reversal and U(l) spin rotational symmetry in 
an integer spin system. The algebraic relations for the 
U(l) x Z2 group are 



TU-gTUe = e. 



(El) 



The corresponding constraints ( 28 ) for symmetry trans- 
formations {W T ,8(j) T } and {W^ = I 2x2l 5(f) U <> = Ot} 
are 



(hx2 - W T )54> T + (J 2 x 2 + W T )9t 



mod 27r, V 




(E2) 



and ( 37|38 ). The gauge inequivalent solutions to these 
constraint equations lead to 




(E3) 



W T = a z , 6</> T = n = 0, 1. (E4) 

\ nir I 

For both n = 0, 1 corresponding symmetry-allowed Higgs 
terms are 



S ldne = J2iez °l I dxdt cos(Z0i) 



a. N = odd integer: Zi classes 

The gauge inequivalent solutions to the above con- 
straint equations are ( 39 I and 



W 9 =I 2x2} Scj> 9 = 



2nk 



kez. 



(E9) 



It's easy to verify that a set of independent Higgs terms 
are J dxdt J^i cos(Z</>i). Hence the variable <f>i can be 
localized at value = without breaking any symme- 
try. So they all correspond to the trivial phase. There 
is only one trivial phase with Zn x Z^ symmetry for 
N = odd. 



b. N — even integer: Minimal set: Z 2 classes 



Solving (|53} and (|E7j) we have W T = a z and W 9 = 

±^2x2- 

(i) For W 9 = I 2x2 the gauge inequivalent solutions are 



n 2 n 



"-'U/J- (E10) 



0< k< N- 1, ui,n 2 = 0,1. 

If tlx = the variable <j>i can be localized at (</>i) — by 
Higgs term — cos <j>\ without breaking the symmetry and 
it corresponds to the trivial phase. If n\ = 1, = k = 
corresponds to the trivial phase again since <f> 2 can be lo- 
calized. Notice that when m = 1 we require (fc, A/2) = 1 
so that transformations (E10 1 form a faithful representa- 
tion of symmetry group G — Zn x Z% , N =even. Let's 



(E5) label a state with the above transformations ( E10 1 as 
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[k, n\, n 2 ] and we have [0, 1, 0] = [k, 0, n 2 ] = ^ ZnxZ t . In 
the following we analyze the group structure formed by 
states [k, n\, n 2 ]. 

Now let's put together a state [k,n\, n 2 ] with edge vari- 
able {4>i,4> 2 } is put together with a state [k', n' x , n 2 ] _1 
with edge variable {(/)[, <p' 2 }, we can choose the following 
independent bosonic variables {k'4>2 — k4> 2 ,k(f>i — k'4>'\} 
and gap all the edge states if (k,k') — 1. The associ- 
ated Higgs terms will preserve the Zn x Z 2 symmetry if 
kn\ — kfn'i = mod 2 and k'n 2 — kn' 2 = mod 2. As a 
result [k, n\, n 2 ]® [k 1 , n\, n'^ 1 — &z n -az t or equivalently 



For (k, k') = 1 : [k,ni,n 2 ] = [k 1 , 



'1: "2J 



if kn\ — k'n'i = mod 2, k'n 2 — kn' 2 = mod 2. 

Therefore we have [2k + l,ni,n 2 ] = [l,ni,n 2 ] by choos- 
ing k' = 1. On the other hand, if k =even we again 
have N/2 =odd since (k, N/2) — 1 for a faithful repre- 
sentation of symmetry group Z^ x Z 2 . We can localize 
bosonic variable {^02 — 2 j T^'i ~ 01 } without break- 
ing any symmetry if we choose k' = 0. Hence we also 
have [2k,ni,n 2 ] — [0, ni,n 2 }. These relations result in 
only three nontrivial SPT phases: [1,1,0], [1,1,1] and 
[0,1, 1]. 

Similarly by putting together a state [fc, l,n 2 ] with 
edge variable {^1,^2} is put together with a state 
[fc', 1,712] with edge variable {4>'\, 4>' 2 }, we can always lo- 
calize bosonic variable cf>i — (/)[ and gap out part of the 
edge. What is left on the edge are described by variables 
{</>i = 0i, 02 = 4 > 2 + 4>' 2 }- They obey Kac-Moody algebra 
(35 1 and transform as a [k + k' , 1, n 2 + n 2 ] state. Hence 
we've shown that 



= [k- 



[fc,l,na]e [k',l,n' 2 ] 
k' mod 2, 1, n 2 + n' 2 



mod 2]. (Ell) 



Since k, n 2 are both Z 2 integers, so clearly all different 4 
states [k, 1, n 2 ] form a Z 2 group. Consequently there are 
3 nontrivial SPT phases labeled by n\ — 1 and [fc, n 2 ] = 
[0,1], [1,0] or [1,1] in fl39t and (|E10|. 

(ii) For W 9 = —I 2X 2 we can always choose a gauge so 
that inequivalent solutions to constraints are 



m,n 2 =0,l. (E12) 

However the above symmetry transformations {W 9 = 
V 



2x2, 







j } do not correspond to a faithful 



representation of Zm x Z 2 group for N =even, unless 
N = 2. And it is not clear how the states with symmetry 
transformations W 9 = — I can be realized in a physical 
bosonic system. Therefore we won't include the states 
with symmetry transformations W 9 = —I in the mini- 
mal set of topological phases with Z^ x Z 2 symmetry, 
as discussed earlier for W 9 = —I phases with Z^ x Z 2 
symmetry. 



In summary, there are Z 2 classes of different bosonic 
non-chiral SRE phase in the presence of symmetry group 
Zn x Z 2 , N =even. When A^=odd there are no nontriv- 
ial SPT phases. 



3. U(l) x Z2 symmetry: Z x Z| classes 

Denoting the group elements of U(l) by Ue, < 9 < 
2tt and generator of Z 2 by g, the group U(l) x Z 2 has 
the following algebraic structure: 



UegU-eg = e. 



(E13) 



and (41 1. The associated constraints for symmetry trans- 



formations are 



{W 9 f = I 2x2 , (J 2x2 + W 9 )60 9 = ri , (E14) 



(I 2x2 + W 9 )5cp 9 + 9(I 2x2 - W 9 )t = r . (E15) 



where we have W Ue = I 2x2 5<p Ue = 9t with t i,t 2 € Z 
and (ti,t 2 ) — 1 for U(l) symmetry. Solving (E14l we 
have W 9 = ±I 2x2 or ±a x . 



a. ZxZ 2 classes with W 9 = ±/2 X 2 

(i) For W 9 — I 2x2 . as guaranteed by theorem (CI I 
we can always transform the "charge vector" t into a 

standard form t = | 



,1) 



q € Z. And the inequiva- 



lent "faithful" symmetry transformations satisfying con- 
straints ( |E14pl5| are 



ni 

,«2, 




ni,n 2 = 0,l, ge(E16) 



Let's label a state with the above transformation rules as 
[q, 711,712]. Similar as earlier discussions for other sym- 
metries, by putting two states [q,ni,n 2 ] and [q', 71^,712] 
together, one can show the following multiplication rule: 



[q,ni,n 2 ] ® W^n'^n^} 
= [q + q', ni + n[ mod 2, n 2 + s 



mod 2] (E17) 



Hence there are Z x (Z 2 ) 2 classes of different phases la- 
beled by integer q and Z2 integers ni,n 2 . The trivial 
phase corresponds to q = n\ = n 2 = 0. 

(ii) For W 9 = —I 2x2 , one can a lways choose a gauge so 

/V 



that 



,0, 



and by solving (El 5) we get t = 



It's easy to verify this corresponds to 



the trivial phase. 
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b. Other solutions to \ElJ\Elh^ with W 9 = ±o\ 



(iii) For W 9 = cr x , the inequivalent solutions to con- 
straint equations are 




, n = 0,l. (E18) 



These two nontrivial SPT phases are labeled as [a x , n] 
with n — 0,1. Their physical realization and group struc- 
ture are not clear. 

(iv) For W 9 = — er x , the inequivalent solutions to con- 
straint equations are 



variables to be localized simultaneously (i.e. they are in- 
dependent bosons) they have to satisfy ( 20 1 and hence 
Y, a A 2 a + B\= 0. This leads to A a = B a = and hence 
no bosonic variable on the edge can be localized without 
breaking the symmetry. Hence whenever we add an ex- 
tra [a z , 0] state into the system, there is one more c = 1 
gapless state on the edge. Hence all the different states 
{[<r x , 0] M © [a, l] N , M, N eZ} form the 1? group, 
(iv) For W 9 = —a z the inequivalent solutions to (99 1 

( o x 



are r\ — and 



where n — 0,1. We label 



1 



-1 



1 



n = 0,1. (E19) 



These two nontrivial SPT phases are labeled as [—a x ,n] 
with n = 0, 1. It's easy to show that 



[a x ,n] 1 = [-a x ,n], n = 0, 1. 



(E20) 



these states by [— o- z ,n] and it's straightforward to show 
that [cr^nj^ 1 = [— cr z ,n]. 

To summarize, with Z 2 symmetry transformation 
W 9 = ±er 2 , there are 1? classes of different fermionic 
non-chiral SRE phases in the presence of Z 2 x Z£ sym- 
metry. It is presently unclear to us if these transforma- 
tion laws can be realized in a physical system of fcrmions. 
We have not found a microscopic realization, hence we 
do not include it in the minimal set of topological phases 
with this symmetry. 



Their physical realization and group structure are not 
clear either as with the discussion in Section IIV D 31 and 
therefore we do not include these phases. 

To summarize, there are Z x (Z2) 2 classes of different 
phases with W T — ±I 2x2 f° r symmetry group U(l) x Z 2 . 

Besides, there are 4 extra possible nontrivial SPT 
phases with W° — ±a x for U(l) x Z 2 symmetry in a 
bosonic non-chiral SRE system. 



Appendix G: Other solutions to (111 I for fermion 
SPT phases with GfjZ{ = Zi x Z'j' symmetry 



W 9 = -I 2 > 



Zo classes 



(ii) If W 9 = —I 2 x2 the gauge inequivalent solutions to 
( 111 I are rj g = and 



Appendix F: Other solutions to (99| for fermion 



SPT phases with Gf/Z( = Z2 symmetry 

Z 2 classes 



1,-9 



±0"z 



r = (f+n)7r 
n,V,VT = 0,1. 




(Gl) 



If Vt = we can destroy the gapless edge states without 
breaking the symmetry, by localizing the bosonic variable 
(j>i — (f)2 at a classical value. If ?; = 0, we can destroy the 
gapless edge states without breaking the symmetry, by 
localizing a differnt bosonic variable <j>i +(^2 at a classical 
value. Hence only 77 = r/x = 1 and n = 0, 1 correspond 
to non-trivial SPT phases, with S(j> T = (tt, 0) t or (0, tt) t 
(n = or 1). Let's label the states with symmetry trans- 
formations (Gil by [rj ) r)T,n]. It's easy to verify that 



(iii) For W 9 = a z the inequivalent solutions to (99) 

(TL7T \ 
where n — 0,1. Since under 

symmetry g we have 0i±02 — ► 0iT02+?i7'", the variables 
cannot be localized without breaking the Z 2 symmetry. 
These two different nontrivial SPT phases are labeled as 
[<7 z ,ti] where n = 0, 1. In the following we identify their 
group structure. 

First notice that when a \a z , 0] state with edge vari- 
ables {0!, 02 } is put together with a [ct^I] -1 state [1, 1, 0] ® [1, T70] = [1, 1, 0] © [1, 1, 1] = e 
with edge variables {^^J,}, its edge cannot be gapped 
without breaking the symmetry, suggesting [a z , 0] 7^ 
[cr z , 1]. Then let's consider TV copies of [a z ,n] states put 
together and their edge variables are {4>^,(t) 2 , 1 < a < 
N} . A generic bosonic variable that can be localized on 
the edge is written as X]^Li(^a^i + B a cj) 2 ), A a , B a g Z 
satisfying J2 a A\ — B\ = due to condition ( fl~9| . Un- 
der Z 2 symmetry generator g this bosonic variable be- 
comes Y^t=i(Aa4'i ~ B a (j>2). In order for the two bosonic 



r f. since 



the edge from two phases put together can be gapped by 
condensing independent bosons {<f>i + (f>' 2 , (f>2 + 4>'\}- Since 
[0,0, n] = [0,1, n] = [1,0, n] = e z 2 xz T xz f we see tnat 
different phases form a Z 2 group. The only nontrivial 
SPT phase is [1, 1, 0] = [1, 1, 1]. The microscopic realiza- 
tion of this particular SPT phase is not clear, (we have 
not found a realization in the coupled wire approach) and 
as discussed previously, we omit it from the minimal set 
of topological phases. 
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bations that tunnel between different layers added to the 
ideal Hamiltonians for multi-layer FQH states (10 1. Since 



there is an energy gap separating the ground states and ex- 
cited states of ideal Hamiltonian (pseudopotentials), any 
weak perturbations that is local will not affect the topo- 
logical properties of the phase. 
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In general if (ti , £2) = t > 1 symmetry transformations 
{W u ° , (f> Ue } will form an unfaithful representation of sym- 
metry group (7(1). Then the symmetry group of ground 
state will be broken from G = U(l) (symmetry group of 
Hamiltonian) down to Zt if the edge is gapped by condens- 
ing bosons in the presence of symmetry-breaking indepen- 
dent Higgs terms. 

These representations however are known to be realized in 
topologically ordered phases with nontrivial ground state 
degeneracy on a torus, such as Z2 spin liquids with trans- 
lational symmetry as discussed in Ref. 1571 



